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The relativistic heavy ion collider (RHIC) was constructed to achieve an asymptotic state of 
nuclear matter in heavy ion collisions, a near-ideal gas of deconfined quarks and gluons denoted 
quark-gluon plasma or QGP. RHIC collisions are indeed very different from the hadronic processes 
observed at the Bevalac and AGS, but high-energy elementary-collision mechanisms are also non- 
hadronic. The two-component model (TCM) combines measured properties of elementary collisions 
with the Glauber eikonal model to provide an alternative asymptotic limit for A-A collisions. RHIC 
data have been interpreted to indicate formation of a strongly- coupled QGP or "perfect liquid." 
In this review I consider the experimental evidence that seems to support such conclusions and 
alternative evidence that may conflict with those conclusions and suggest different interpretations. 

PACS numbers: 12.38. Qk, 13.87.Fh, 25.75.Ag, 25.75.Bh, 25.75.Ld, 25.75.Nq 



I. INTRODUCTION 

The experimental program addressing high-energy 
heavy ion (A-A) collisions initiated at the Berkeley Be- 
valac and later extended to the Brookhaven alternating 
gradient synchrotron (AGS) reported collective motion 
(flow) of nucleon clusters, nucleons and hadronic res- 
onances for center-of-mass collision energies below 10 
GeV [Tj. Arguing by analogy with the thermodynam- 
ics of ordinary matter and QCD asymptotic freedom it 
was conjectured that increased beam energy might lead 
to energy densities and temperatures sufficiently high 
to "melt" colorless hadrons to produce a thermalized 
medium of deconfined weakly-coupled colored quarks and 
gluons referred to as a quark-gluon plasma (QGP) [2]. 
Detection of collective flow in the produced medium 
might support the inference of thermalization required 
to conclude that a new form of QCD matter had been 
achieved. The Brookhaven relativistic heavy ion collider 
(RHIC) was proposed to achieve such energy densities 
and to produce a QGP. 

A different context (HEP) emerges from observations 
of high-energy elementary hadronic collisions and their 
QCD description. Two mechanisms appear to domi- 
nate high-energy p-p collisions: (a) projectile dissocia- 
tion to low-transverse-momentum p t hadrons extending 
over a broad longitudinal rapidity y z interval according 
to a parton distribution function (PDF) and (b) infre- 
quent large-angle scattering of some partons according to 
perturbative QCD (pQCD) followed by nonperturbative 
fragmentation to collimated jets of hadrons with a broad 
distribution of momenta, possibly extending to large as 
well as small values. That scenario is the two-component 
(soft+hard) model (TCM) of hadronic collisions. 

Prior to acquisition of RHIC data the possibility ex- 
isted that A-A collisions might be just linear superpo- 
sitions of N-N collisions following the TCM scenario. 
The null hypothesis for hadron production would then 
be fragmentation, either longitudinal (soft) or transverse 
(hard). The HEP limiting case balances the QGP alter- 
native in which most hadrons emerge from a thermalized 



quark-gluon medium by a little-understood process. The 
reality of high-energy A-A collisions might lie somewhere 
between, depending on imposed conditions, and analysis 
methods should anticipate the full range of possibilities. 

The RHIC theoretical and experimental programs have 
evolved to emphasize the context extrapolated from lower 
energies. A QGP is formed by early thermalization of a 
gluon condensate (color glass condensate or CGC) also 
called a glasma [3] . Hadron production is dominated by 
freezeout from the resulting flowing bulk medium com- 
posed of quarks and gluons. Jet manifestations are con- 
fined to a small fraction of all hadrons, those emerging 
as fragments from high-energy scattered partons. 

The hadron transverse momentum spectrum is accord- 
ingly divided into three intervals. Low p t 0-2 GeV/c 
should include manifestations of soft processes within the 
bulk medium including radial and elliptic flows. Statis- 
tical models of hadron species abundances are applica- 
ble to the thermalized system. Fluctuations of conserved 
quantities may reveal proximity to a phase boundary sep- 
arating hadronic and partonic phases. Hydrodynamic de- 
scriptions of flows may reveal a QCD equation of state. 

Intermediate p t 2-5 GeV/c is assumed to represent a 
transition interval between bulk medium freezeout and 
pQCD jet production. The theoretical description is 
uncertain and may be driven by data phenomenology. 
Quark recombination or coalescence models applied to 
jet-medium interactions are proposed. 

High p t above 5 GeV/c should be dominated by 
hadrons from jets. A pQCD description of parton scatter- 
ing and fragmentation to jets is applicable. Evidence for 
modified jet production (jet quenching) resulting from 
parton energy loss and other interactions in the dense 
colored medium is sought. Isolation of jets from a large 
combinatoric background presents technical challenges. 

After three years of data acquisition and analysis the 
theoretical consensus was developed that not a weakly- 
coupled but a strongly- coupled QGP or sQGP is formed in 
central 200 GeV Au-Au collisions. The material appar- 
ently exhibits a very small viscosity and little dissipation, 
suggesting the term "perfect liquid." 
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The experimental and theoretical output from RHIC 
to date is quite large. The scope of this review is re- 
stricted as follows. Topics dealing directly with claims 
for formation of a QGP or "perfect liquid" are preferred. 
Results from three of the four experiments, those em- 
phasizing the central rapidity region (collision center of 
momentum) are preferred. A small sample of figures and 
papers for each topic is presented, those seeming to illus- 
trate major results for the topic. Conventional analysis 
according to assumptions outlined above is compare with 
alternative analysis that may be more compatible with 
the HEP context. 

This paper is organized as follows: Sec.[H|reviews con- 
ventional and alternative analysis methods, Sec. |III| re- 
views results from the first three years of RHIC oper- 
ations, Sec. |IV| summarizes white papers from the four 
RHIC collaborations, Sec. [V] reviews theoretical argu- 
ments leading to claims for formation of a strongly- 
coupled QGP or "perfect liquid." Sec. |Vl] reviews sub- 
sequent progress on spectrum analysis, Sec. |VII| su m- 
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marizes results for elliptic flow analysis, Sec. 
views jet studies and Sec. |IX| reviews fluctuation analysis. 
Sec. [X]pr esents results from alternative analysis methods. 
Sees. |Xl] and |XIl| present discussion and summary. 



II. ANALYSIS METHODS 

Many analysis methods developed for RHIC data prior 
to start-up assumed an established theoretical context. 
The A-A final state should emerge from an expanding 
bulk medium, the result of "freezeout" or decoupling of 
that medium to free-streaming particles. The proper- 
ties of the medium, whether a hadronic or quark-gluon 
scenario better described it, were of primary interest. 
Access to medium properties should be through hadron 
yields and spectrum analysis, with the so-called blast- 
wave model preferred for spectra. Bulk properties should 
be accessed primarily at low p t < 2 GeV/c. 

Some properties of the medium might also be probed 
by high-energy scattered partons fragmenting to jets. 
Would some jets be absorbed by the medium, would 
some jets survive with modification by the medium? How 
would the medium in turn be modified by jets? Jet infor- 
mation and theoretical definition should be best accessed 
at high p t > 5 GeV/c where a pQCD description of jet 
production and jet fragmentation (DGLAP) is valid, and 
jet fragments dominate the soft component of spectra. 

That QGP scenario is an alternative to the high-energy 
physics description of A-A collisions in terms of "direct 
reactions" dominated by projectile-nucleon dissociation 
and in-vacuum parton scattering to jets with fragment 
distributions extending down to small hadron momenta. 
In the latter case a different combination of analysis 
methods might be more appropriate. 



A. Conventional methods 

A-A collision centrality is measured by charge mul- 
tiplicity density dn c h/dr) or participant- nucleon number 
Np ar t which emphasizes the upper half of the A-A total 
cross section and tends to obscure the lower half. Spectra 
are presented on transverse momentum p t which tends to 
emphasize structure at higher p t for a given p t acceptance 
interval. Spectra from more-central Au-Au collisions are 
compared directly to p-p or d-Au collisions. There is no 
reference model for transparent A-A collisions (i.e., for no 
partonic or hadronic secondary interactions). Jet mod- 
ification is measured by spectrum ratios Raa or Rcp 
which include soft as well as hard spectrum components. 
Sensitivity to jet modifications (hard component) is then 
confined to p t > 4 GeV/c. 

Structure (anisotropy) in ID azimuth correlations is 
measured primarily by V2, the square root of a per-pair 
correlation measure (defined below) , interpreted to repre- 
sent elliptic flow. Other correlation structure (e.g., rep- 
resenting jets) is isolated by subtracting a background 
based on numerical V2 measurements. Jet-related ID 
and 2D angular correlations are further conditioned by 
trigger-associated pt cuts based on assumptions about the 
relative contribution of jets to spectra and correlations. 
Jets may be assumed to retain the same form as in p-p 
collisions. Any deviations from the p-p jet form are at- 
tributed to nonjet contributions (e.g., a "ridge"). Fluc- 
tuations are measured by per-pair statistical quantities 
that include a trivial but dominating l/n c h or l/N part 
trend. 



B. Alternative methods 

A-A collision centrality is measured by fractional cross 
section a~/o~a, fractional impact parameter b/bo or mean 
participant-nucleon path length v = 2Nt> in /N part . Those 
centrality measures provide balanced comparisons of pe- 
ripheral and central collision systematics over the full 
centrality range down to N-N (nucleon-nucleon) colli- 
sions. Spectra are plotted on transverse rapidity y t = 
ln[(mt J rPt)/rn} 1 \ (rn/j is the hadron mass, with default 
TOjr) which provides more-balanced visual comparison of 
low-p t vs high-pi structure. Jet modifications are mea- 
sured by differential ratio taa which compares only the 
hard components of spectra, providing significant access 
to jet structure and modifications down to 0.5 GeV/c [4]. 

Correlations are measured with per-particle measure 
Ap/^/p re f which is by definition independent of system 
size under combination of un correlated subsystems [5]. 
Four charge combinations are studied: LS (like-sign), US 
(unlike-sign) , CI = LS + US (charge- independent), CD 
= LS — US (charge-dependent). Significant features in 
correlations are modeled by simple functional forms not 
motivated by a priori physical models [6] . 

Initial studies focus on minimum-bias correlations (no 
Pt cuts). Thus, 100% of the jet structure for all centrali- 
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ties is considered. A Glauber linear superposition (GLS) 
reference (transparent A-A collisions) is defined based on 
p-p measurements and the Glauber model of A-A colli- 
sions. Significant deviations from GLS are emphasized. 

A nonjet quadrupole component (cos[2((^i — <fe)]) is iso- 
lated from jet structure by 2D model fits. A quadrupole 
spectrum inferred from its p t dependence can be com- 
pared directly with hydro boost predictions. Fluctua- 
tions are also measured with a per-particle quantity. The 
scale (bin-size) dependence of fluctuations so measured 
can be inverted to infer underlying angular correlations 
that are directly and simply comparable with theory. 



III. THE FIRST THREE YEARS 

The main goal of the RHIC in 2000 remained detection 
and study of the QGP as a near-ideal gas of quarks and 
gluons by observing certain "signals" in heavy ion colli- 
sion data. A variety of analysis methods was prepared 
to search for those signals. Analysis methods were tested 
on large volumes of Monte Carlo data prior to first RHIC 
data in August, 2000. In this section we review some of 
the more prominent early results that seemed to support 
claims for formation of a QGP in Au-Au collisions. 



A. Hadron production 

Total charge multiplicity n c h for most-central A-A col- 
lisions, the event-number minimum-bias distribution on 
rich and the charge multiplicity variation with centrality 
provide initial indications of new physics that might be 
accessed by a higher collision energy. 

Figure [I] (left) shows charge multiplicity per 
participant-nucleon pair vs N part near 77 = for 200 and 
62 GeV Au-Au collisions compared to theoretical predic- 
tions [7]. The data span the most-central 40% of the 
Au-Au total cross section and increase relative to a p-p 
value near 2.4 to a value near 3.8 for central collisions. 
The theoretical curves represent a two-component (mini- 
jet) model [5] and a saturation-scale or CGC model 0. 

Figure [l] (right) shows the ratio of 200 to 62 GeV yield 
data compared to HIJING [TO] results representing the 
two-component model. The conventional HEP descrip- 
tion of A-A collisions assuming linear superposition of 
N-N collisions includes minijet production in a TCM. 
The HIJING Monte Carlo is elected to represent that 
model [IT. The different trends in Fig. [I] (right) are 
then interpreted to reject the TCM in favor of a color 
glass condensate (CGC) model [3]. Such yield trends are 
considered essential to establish initial conditions for hy- 
drodynamic theory calculations and to test the CGC hy- 
pothesis. The apparent deviation of the hadron produc- 
tion data from HIJING is characterized as a "too slow" 
increase of the hadron yield with increasing centrality 
and collision energy [TOJ [TO] . 
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FIG. 1: Left: Charge multiplicity density per participant 
nucleon vs N part for 130 and 200 GeV Au-Au collisions 
(points) compared to theoretical predictions (curves) [?]■ The 
lower panel indicates that the energy ratio does not vary 
significantly with centrality. Right: The energy ratio for 
data (points) compared to that predicted from the HIJING 
Monte Carlo (curve) based on minijet production in a two- 
component model, showing apparent disagreement. 



B. Low-pt spectra and radial flow 

In the \ow-pt interval 0-2 GeV/c particle production 
is assumed to be dominated by freezeout (chemical and 
kinetic) from a flowing bulk medium, possibly a thermal- 
ized QGP. The statistical model should predict hadron 
species abundances corresponding to a chemical freeze- 
out temperature T c ^ em and relevant chemical potentials 
fi. Radial flow should be manifested as a deviation of 
p t or m t spectra from a Maxwell-Boltzmann reference 
function. So-called "blast-wave" fits are used to model 
spectra with parameters Tki n , the kinetic freezeout tem- 
perature and (j3t), the mean radial speed. 
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FIG. 2: Bulk-medium parameters inferred from statistical 
model and blast-wave model fits to yields and pt spectra plot- 
ted vs total charge density for identified hadrons from 200 
GeV Au-Au collisions [14] . Left: Fitted chemical T c h and ki- 
netic Tkin freezeout temperatures. Right: Fitted mean trans- 
verse speed {Pt)- 

Figure [2] (left) shows chemical T c h and kinetic T^n 
freezeout temperatures (open points) vs Au-Au centrality 
(charge multiplicity) inferred from blast-wave fits to PID 
spectra for pions, kaons and protons below m t — = 0.6 
GeV/c 2 Q3]. The fit results are typical: T ch remains near 
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150 MeV approximately independent of centrality. Tki n 
falls from a higher value near 130 MeV for peripheral col- 
lisions to a lower value near 90 MeV for central collisions. 

Figure [2] (right) shows mean radial speed (/3t) inferred 
from the same fits. The data increase from 0.25 for p-p 
collisions to about 0.6 for central Au-Au collisions. The 
combination of temperature and radial speed trends is 
interpreted to indicate isentropic expansion of a thermal- 
ized medium. The hadron species ratios are determined 
earlier, and the hadronic system then continues to do 
thermodynamic work until kinetic freezeout at a later 
time. 



C. Elliptic flow 



deviate from expectations for ideal hydrodynamics. The 
jet-like correlation structure may represent the first di- 
rect evidence for jets at RHIC. 

Figure [4] (left) compares ratio v 2 /e from the AGS, SPS 
and RHIC vs hadron density that together seem to con- 
firm a conjecture in Ref. |18j : The v 2 /t ratio (possi- 
bly measuring transport of initial-state nucleon spatial 
anisotropy to final-state hadron momentum anisotropy) 
should increase from near zero for peripheral collisions 
and low hadron densities (little secondary scattering) 
to a large value (possible saturation) for central colli- 
sions and a possible QGP [12] ■ v 2 /e values for most- 
central RHIC collisions appear to achieve a predicted 
ideal- hydrodynamic limit (band). 



"Elliptic flow" is the interpretation of the correlation 
structure represented by measure v 2 = (cos[2(0 — & r )]) 
inferred from two-particle correlations on azimuth rela- 
tive to estimated reaction-plane angle The magni- 
tude of elliptic flow at RHIC in comparison to lower- 
energy data from the SPS has been of major interest. 
The v 2 values at the SPS fell below hydrodynamic pre- 
dictions. Would (ideal) hydro apply to RHIC collisions? 
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3: Elliptic flow V2 data from 130 GeV Au-Au colli- 
sions [T5] . Left: p t -integral V2 vs centrality measured by total 
charge multiplicity relative to the maximum observed value. 
Right: Minimum-bias V2 vs transverse momentum p t . 

Figure [3] shows the first v 2 measurements on centrality 
and pt (and the first paper on any topic) from RHIC [15] . 
The surprising result was the near equivalence of v 2 data 
(solid points) and (ideal) hydro theory (open boxes). The 
data also increase strongly with hadron p t (right panel). 
The V2{pt) results are minimum-bias (averaged over cen- 
trality) . In Ref. |16j v 2 data are inferred from Fourier se- 
ries fits to two-particle azimuth correlations and plotted 
vs fractional cross section for several pt bins. The data 
are interpreted to indicate transformation of initial-state 
A-A overlap geometry to final-state momentum asymme- 
try. Deviations oiv 2 /e scaling at larger p t are interpreted 
to indicate a possible change in the EoS. 

A similar analysis in Ref. [17] reported v 2 (pt) results 
and apparent jet structure derived from azimuth corre- 
lations. Possible saturation of v 2 above 3 GeV/c would 
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FIG. 4: Left: Ratio V2/e vs particle density (1 / S)dn c h / 'dr\ for 
several collision energies, where S is the A-A overlap area [19] . 
The data are compared with an ideal-hydro prediction (band) . 
Right: Distribution of V2 on 77 showing a sharply-peaked dis- 
tribution [201 . 



Figure [4] (right) shows v 2 variation over a large T) in- 
terval (averaged over centrality) (50] . Hydro theory has 
difficulty with such a trend (strong variation on 77) [21] . 
And v 2 (what is actually inferred from two-particle cor- 
relations) would be even more sharply peaked. 

A critical test of the hydro interpretation of v 2 data is 
PID data for identified hadrons. If hadrons emerge from 
a radially-expanding bulk medium then "mass ordering" 
(hadron mass increase from left to right on p t ) should 
be observed below 2 GeV/c. The same mass trend is 
expected for mt spectra as a manifestation of radial flow. 

Figure [5] (left) shows minimum-bias PID v 2 data from 
Ref. [53]. The requisite mass ordering below 2 GeV/c 
is observed, seeming to confirm a hydro interpretation. 
Above 2 GeV/c the data (points) deviate strongly from 
ideal hydro theory (curves), possibly indicating satu- 
ration. Fig. [5] (right) shows similar spectra for sev- 
eral Au-Au centrality bins and for (anti)particles at 
(top)bottom. Although the amplitude trend increases 
with decreasing centrality (as for p t -integral data) the 
trends on p t are very similar over a large centrality and 
Pt interval. 
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FIG. 5: Left: V2 vs p t for identified hadrons from minimum- 
bias 200 GeV Au-Au collisions (points) compared to ideal 
hydro predictions (curves) [23]. Possible saturation of data 
above 2 GeV/c is indicated. Right: PID i>2 (pt) data for sev- 
eral Au-Au centrality intervals showing similar structure [24] . 



D. High-pt spectra and jet correlations 

If a color-deconfined medium is formed in A-A col- 
lisions scattered high-energy partons, as QCD colored 
objects, should interact with the medium and lose en- 
ergy. The result is termed jet quenching: most jets are 
modified and some jets may be entirely absorbed by the 
medium. Jets should be manifested in single-particle pt 
spectra (by assumption only at higher p t ) and in two- 
particle angular correlations. The reference for spectrum 
analysis is p t spectra from p-p collisions. The spectrum 
ratio R AA = (l/N bin )p AA (p t )/ p pp (p t ) compares jet man- 
ifestations in A-A collisions to in vacuum jets in p-p col- 
lisions. In case of linear superposition of N-N collisions 
within A-A collisions R AA should be 1 over a pt interval 
where dijet production dominates spectra. Values sub- 
stantially less than 1 would indicate jet quenching. 

Figure [6] (left) shows data from Ref . [35] (early results 
for R AA ) that reveal values of R AA substantially less 
than 1 above 2 GeV/c, seeming to confirm substantial jet 
quenching in more-central Au-Au collisions. The RHIC 
results are compared with Pb-Pb results from the SPS 
at 17 GeV (solid curves) that exceed 1 as expected for 
initial-state effects (Cronin effect) and no jet quenching. 

Figure [6] (right) shows a comparison between Au-Au 
(lower points, central) and d-Au (upper points) showing 
that partons passing through cold nuclear matter (the 
latter case) do not manifest jet quenching [26] . It is con- 
cluded that jet quenching is only observed for the hot and 
dense medium created in more-central Au-Au collisions. 

Figure [7] (left) from Ref. [27] shows R AA spectrum 
ratios vs Au-Au centrality from central (upper left) to 
peripheral (lower right). Jet quenching proceeds from 
strong (R AA » 0.2) to negligible {R AA ~ 1) over that 
centrality range, apparently illustrating the centrality 
evolution of QCD medium formation. Similar results 
were reported in Ref. [2"5] . 

Figure [7] (right) from Ref. 26j shows dihadron correla- 
tions from p-p, d-Au and Au-Au collisions. A narrow 
same-side peak representing mirajet correlations from 




FIG. 6: Left: Spectrum ratio Raa vs pt for 130 GeV 
Au-Au collisions (points) showing suppression (jet quenching) 
above 2 GeV/c in marked contrast to substantial enhance- 
ment (Cronin effect) for Pb-Pb collisions at 17 GeV (solid 
curves). [25] Right: Comparison between strong suppression 
in central 200 GeV Au-Au collisions (lower points) and no 
suppression in 200 GeV d-Au collisions (upper points), con- 
firming no jet suppression in cold nuclear matter [26] 
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FIG. 7: Left: Centrality dependence of R A A(pt) from 200 
GeV Au-Au collisions [27] ■ Strong suppression for central 
collisions reduces to negligible suppression for peripheral col- 
lisions. Right: Jet-related dihadron azimuth correlations 
demonstrating disappearance of the away-side jet (A(j> — n) 
in central 200 GeV Au-Au collisions (blue points) compared 
its presence for p-p collisions (histogram) 26 . The away-side 
peak is also present in d-Au collisions (red and green points) . 



"triggered" jets is expected and observed at A<fi = for 
all cases. The critical question is whether an internet 
away-side peak representing recoil-partner jets appears 
at A(j) — 7T. For p-p, d-Au and peripheral Au-Au colli- 
sions an away-side peak is observed representing the re- 
coil partner of a dijet pair. For central Au-Au collisions 
the peak representing the away-side jet has disappeared. 

Such high-p t spectrum and correlation results are in- 
terpreted to demonstrate formation of a dense colored 
medium in which some fraction of jets is fully absorbed 
(opaque medium). Triggered jets are believed to be bi- 
ased to emerge only near the medium surface (surface 
bias), and the recoil partners are then biased to pass 
through most of the dense medium and lose much or all 
of their energy in more-central Au-Au collisions. The evi- 
dence from R AA and dihadron correlation data in several 
collision systems is seen as consistent with that scenario. 
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There is apparently 80% reduction of all jet fragments 
in central Au-Au collisions, and triggered jets have no 
detectable recoil partner given the imposed p t conditions. 



E. Minimum-bias 2D angular correlations 

Although ID dihadron analysis on azimuth imposes 
specific trigger-associated p t cuts considered appropriate 
for high-p t jet analysis it is possible to study angular cor- 
relations on both pseudorapidity and azimuth with no p t 
conditions, resulting in minimum-bias or p t -integral 2D 
angular correlations. There are no a priori expectations 
for correlation structure. For a thermalized system there 
might be no significant jet correlations when averaged 
over all combinatoric pairs. The per-particle measure 
N(r — 1) in the following plots is related to the later 
per-particle measure Ap/ \fp re f [5] by 
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i/NIf 2irAi] 



(r-1) 



(1) 



where N/2ttAi] w po(b) is the average single-particle 2D 
angular density near midrapidity, as noted in [5|. The 
updated measure is a density approximately independent 
of detector acceptance near mid rapidity. 

Figure [8] shows CI correlations for peripheral (left) and 
central (right) 130 GeV Au-Au collisions [6]. Angular 
difference variables (e.g., t^a = f]i — t\i and similar for 
azimuth) are appropriate for such untriggered combina- 
toric correlations where no trigger-particle reference an- 
gles apply. Those minimum-bias 2D angular correlations 
include three prominent structures: (a) a SS 2D peak 
modeled by a 2D Gaussian, (b) a broad AS ID peak on 
azimuth modeled by dipole cos((/>a — 7r) and (c) a nonjet 
azimuth quadrupole modeled by cos(20a)- In Fig. [8 fitted 
features (b) and (c) have been subtracted to emphasize 
evolution of the SS 2D peak with centrality. 




FIG. 8: Minimum-bias (p t -integral) CI (charge-independent) 
2D angular correlations on angle differences (77a><^a) from 
130 GeV Au-Au collisions [6]. A fitted broad away-side ID 
peak at <^a = 7r and nonjet quadrupole cos(20a) have been 
subtracted. The remaining same-side 2D peak is strongly 
elongated proceeding from peripheral (left, consistent with 
p-p data [3T]) to central (right) Au-Au collisions. 



Correlation features (a) and (b) are expected for dijet 
pairs: intrajet correlations (a) + interjet correlations (b). 
But the correlations in Fig. [8] were formed from particles 
with pt < 2 GeV/c where hydro (elliptic and radial flows) 
is expected to describe spectra and correlations. The jet- 
like result was quite surprising. Feature (c) is nominally 
consistent with elliptic flow expectations, and its ampli- 
tude serves as an alternative measure of t>2 as inferred 



from model fits to 2D angular correlations (Sec. X D ) 



The narrow peaks at the origin represent conversion elec- 
tron pairs and Bose-Einstein (quantum) correlations. 

Along with the surprise of large-amplitude jet-like an- 
gular correlations below 2 GeV/c (Sec. IV A I came evolu- 



tion of the SS 2D peak shape from strongly elongated on 
<j> (p-p and peripheral Au-Au) [301 EI] to strongly elon- 
gated on 77 (more-central Au-Au) [6] . The 77 elongation 
(so-called same-side "ridge" ) has become a central prob- 
lem for theory and experiment in later RHIC and LHC 
analysis. 

The plots in Fig. [8] result from subtracting a fitted 
nonjet quadrupole and a large-amplitude AS ID peak 
fully consistent with back-to-back jet correlations rela- 
tive to the SS 2D peak. Given that interpretation the 
"away-side" peak representing jet-jet correlations does 
not disappear in central Au-Au collisions but is signif- 
icantly modified (mean p t is reduced). A similar AS 
dipole is expected for combinatoric hadronic (not dijet) 
momentum conservation, but the amplitude and central- 
ity dependence for that feature are very different. If the 
combination (a) + (b) is indeed (mini)jets then argu- 
ments against a two-component minijet model based on 
HIJING simulations [TT] can be questioned (Sec. XI B). 

Figure [9] shows complementary CD angular correla- 
tions from the same 130 GeV Au-Au collision data for pe- 
ripheral (left) and central (right) collisions. Negative val- 
ues reflect "canonical suppression" of net charge within 
small angle differences. A large-amplitude 2D peak at 
the angular origin may be interpreted as representing lo- 
cal charge conservation during parton fragmentation to 
(mini)jets (local charge conservation within a fragmenta- 
tion cascade) 32, 33]. For peripheral collisions (left) an 
additional smaller-amplitude broad ID peak on ?7a can 
be interpreted in terms of projectile-nucleon fragmenta- 
tion to charge-neutral hadron (mainly pion) pairs. For 
central collisions a small-amplitude structure narrow on 
</>A and extended on 77A corresponds to r\ elongation of 
the SS 2D peak ("ridge"). 

The charge balance function (BF) was proposed to de- 
tect "delayed hadronization" from a partonic medium by 
broadening of CD angular correlations on r\. The predic- 
tion is based on arguments relating to longitudinal ex- 
pansion of the QCD medium [34] . It was shown that the 
BF is a ID projection of CD 2D correlations in Fig. [9] 
onto ID ?7a |35j . The narrow structure on azimuth is 
lost in such projections. The possible dominance of jet 
contributions to CD structure is then overlooked. These 
CD 2D angular correlations are consistent with the TCM 
for hadron production: longitudinal projectile fragmenta- 



FIG. 9: Left: Minimum-bias (pt-integral) CD (charge- 
dependent) 2D angular correlations from 130 GeV Au-Au 
collisions [33]. A narrow 2D exponential peak at the origin 
is common to all centralities. A broader ID peak on rjA is 
observed in more-peripheral collisions (left). An extension 
narrow on 0a is observed in more-central collisions (right) 
that may be associated with the r\ elongation in Fig. [5] 



FIG. 10: Left: Mean-pt fluctuations measured by r.m.s. 
quantity A(j pt:n vs centrality measured by n c h from 130 GeV 
Au-Au collisions [36]. Upper points represent uncorrected 
CI (pt) fluctuations, lower points represent CD (related to 
net-charge) fluctuations. The solid curves are corrected for 
detector inefficiencies. Right: Mean-p t fluctuations measured 
by alternative r.m.s. measure F vt vs centrality measured by 
Npart showing a similar non-monotonic trend [38] . 



tion and transverse parton fragmentation, both with local 
charge conservation. There is no significant CD structure 
corresponding to the nonjet azimuth quadrupole promi- 
nent in CI correlations. 



F. p t and number fluctuations 

If RHIC A- A collisions achieved sufficient energy densi- 
ties and temperatures to cross the QCD phase boundary 
and form a QGP event-wise critical fluctuations in some 
event properties might be detectable based on analogies 
with macroscopic critical phenomena. A candidate event 
property is event-wise mean p t denoted (pt) interpreted 
as a proxy for kinetic temperature T^ in . The question is 
posed: Do RHIC collisions exhibit critical temperature 
fluctuations for some collision conditions? 

(left) shows fluctuation measure Acr Pt:n , a 



Figure 
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per-particle r.m.s. measure of mean-p 4 fluctuations, vs 
centrality measured by relative charge multiplicity N/Nq, 
where Nq estimates the mean multiplicity corresponding 
to b = [37], for 130 GeV Au-Au collisions [SB]. Fluc- 
tuations are measured at the "acceptance scale," the an- 
gular acceptance of the STAR TPC. The positive points 
correspond to CI fluctuations, the negative to CD (net- 
charge) fluctuations. The centrality variation is smooth 
from very peripheral to very central collisions. There 
is no discontinuity suggesting a phase boundary. The 
fluctuation amplitude is large, prompting the question — 
what mechanism can produce such large (p t ) fluctuations 
over all Au-Au centralities? 

(right) shows mean-p t fluctuations mea- 



Figure 10 



sured by an alternative statistic F p 



A relative- 



variance measure uj pt — a^ p ^/pf is evaluated for sin- 
gle events and for a mixed-event reference. Ratio r — 
u Pt , single /u P t, mixed is defined and F Pt = r- 1. The same 
general trend of increase then decrease with centrality 
is followed, but the numerical values are quite different. 



The RHIC and SPS communities were confronted with a 
variety of fluctuation measures with very different prop- 
erties and interpretations [39) . The lack of conventions 
for fluctuation measures, the dependence on specific de- 
tector geometries and imposed kinematic conditions, re- 
mained an unresolved problem. More-recent fluctuation 
measures and interpretations are discussed in Sec. |XE[ 
In general, no sharp variations in fluctuation trends that 
might correspond (by analogy) to a macroscopic phase 
transition were observed in early RHIC fluctuation data. 



IV. THE RHIC WHITE PAPERS 

During 2004 the four RHIC experiments reviewed the 
analysis performed and data obtained during the first 
three years of RHIC operation. The review results were 
presented as white papers [4"UH4"3"] . The central question: 
had convincing evidence emerged to support discovery 
of a QGP? If not, what evidence was missing and what 
might be required to obtain it. If so, what further studies 
were required for its complete description. In this section 
I present brief highlights from each of the RHIC white 
papers. 



A. STAR 

The QGP is first defined: "...a (locally) thermally equi- 
librated state of matter in which quarks and gluons are 
deconfined from hadrons, so that color degrees of freedom 
become manifest over nuclear, rather than merely nucle- 
onic, volumes." The expression "strongly-interacting" 
QGP is said not to differentiate from that expected in 
the temperature regime just above T c accessed by RHIC 
collisions, a highly-correlated quark-gluon system. 

Standards of scientific proof are identified: necessity - 
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do data require a QGP explanation with no other alterna- 
tives; consistency - are specialized theories applied to dif- 
ferent data aspects mutually compatible; endurance - are 
incompatibilities with some data indicative of inevitable 
small adjustments or a failure of the basic paradigm? 

Experimental results and relevant theories are re- 
viewed for hadron production, hadron species abun- 
dances, spectrum structure (flows and jets), elliptic flow 
V2 systematics, fluctuations and jet-related correlations. 

Resulting argument in favor of QGP includes "...the 
unprecedented success of hydrodynamics in providing 
a reasonable quantitative account for collective flow at 
RHIC, and of the statistical model in reproducing hadron 
yields through the strange sector, together argue for an 
early approach toward thermalization spanning the u, d 
and s [quark] sectors." 

Also, "The fitted parameters of the statistical model 
analyses, combined with inferences from the produced 
transverse energy per unit rapidity, suggest attainment of 
temperatures and energy densities at least comparable to 
the critical values [e.g., T c ] for QGP formation in LQCD 
calculations of bulk, static strongly interacting matter," 
implying that colored quarks and gluons are the relevant 
dynamical degrees of freedom. 

However, the paper cautions: "On the other hand, 
measurements of angle difference distributions for soft 
hadron pairs [e.g., Fig. H of this paper] reveal that some 
(admittedly heavily modified) remnants of jetlike dynam- 
ical correlations survive the thermalization process, and 
indicate its incompleteness." 

Among several prominent goals for future work are: 

(a) establish that jet quenching is an indicator of par- 
ton, and not hadron, energy loss and (b) extend RHIC 
Au-Au measurements down toward SPS measurements 
in energy, to test quantitative predictions of the energy 
dependence. At present (a) has been accomplished and 

(b) is in progress as the RHIC beam energy scan (BES). 



B. PHENIX 

Based on a summary of experimental results it is con- 
cluded that "...hydrodynamical calculations which repro- 
duce the magnitude of elliptic flow observed at RHIC re- 
quire local thermalization to occur very quickly, typically 
by 1 fm/c or earlier.... If the system does reach local equi- 
librium on this time scale then the energy density of the 
first thermalized state would be in excess of 5 GeV/fm 3 , 
well above the amount required to create the QGP." 

Concerning minijets and the TCM or HEP context: 
"Pre-RHIC expectations that E t and charged particle 
production would be dominated by factorized pQCD pro- 
cesses [jets] were contradicted by data, which showed only 
very modest increases with centrality and beam energy. 
A new class of models featuring initial-state gluon satura- 
tion compares well with RHIC multiplicity and E t data, 
and are also consistent with our Bjorken style arguments 



for estimating energy densities at early times." The ob- 



servation and its implications are discussed in Sec. XI B 



"The observed suppression of high-p t particle produc- 
tion at RHIC... provides direct evidence that Au+Au 
collisions at RHIC have produced matter at extreme den- 
sities, greater than ten times the energy density of nor- 
mal nuclear matter and the highest energy densities ever 
achieved in the laboratory." 

The paper concludes that "...there is compelling exper- 
imental evidence that heavy-ion collisions at RHIC pro- 
duce a state of matter characterized by very high energy 
densities..., strong collective flow, and early thermaliza- 
tion. ...this matter modifies jet fragmentation and has 
opacity that is too large to be explained by any known 
hadronic processes. [...] The most economical description 
is in terms of the underlying quark and gluon degrees of 
freedom." 



C. PHOBOS 

Emphasis is placed on particle production and the 
large rj acceptance of the PHOBOS detector. Inferred 
scaling behaviors, factorizations (e.g., centrality and en- 
ergy) and manifestations of the initial-state A-A geome- 
try in the data are featured. Descriptions of the data in 
terms of hadronic degrees of freedom are ruled out. 

"...the data clearly demonstrate that proportionality 
to the number of participating nucleons, N part , is a key 
concept which describes much of the phenomenology [...] 
Further, the total particle yields per participant from dif- 
ferent systems are close to identical when compared at 
the same available energy.... [...] ...and many charac- 
teristics of the produced particles factorize to a surpris- 
ing degree into separate dependences on centrality and 
beam energy. [...] All of these observations point to the 
importance of the geometry of the initial state and the 
very early evolution of the colliding system in determin- 
ing many of the properties of the final observables." 

The paper concludes "...these simple scaling features 
will constitute an integral component or essential test of 
models which attempt to describe the heavy ion collision 
data at ultrarelativistic energies. These unifying features 
may, in fact, provide some of the most significant inputs 
to aid the understanding of QCD matter in the region 
of the phase diagram where a very high energy density 
medium is created." 



D. BRAHMS 

Concerning evidence for the QGP in data it is pointed 
out that "...experimental signatures can be roughly 
grouped into two classes: 1) evidence for bulk properties 
consistent with QGP formation, e.g. large energy den- 
sity, entropy growth, plateau behavior of the thermody- 
namic variables, unusual expansion and lifetime proper- 
ties of the system, presence of thermodynamic equilibra- 
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tion, fluctuations of particle number or charge balance 
etc, and 2) evidence for modifications of specific proper- 
ties of particles thought to arise from their interactions 
with a QGP, e.g. the modification of widths and masses 
of resonances, modification of particle production proba- 
bilities due to color screening (e.g. J/^B suppression) and 
modification of parton properties due to interaction with 
other partons in a dense medium (e.g. jet quenching), 
etc." 

The paper asks "...if there are any specific features that 
may falsify the conclusion [of QGP formation]? To our 
knowledge no tests have been proposed that may allow 
falsification of cither a partonic scenario or a hadronic 
scenario, but it would be important if any such exclusive 
tests were to be formulated." 

Based on a summary of available data the paper con- 
cludes "The overall scenario is therefore consistent with 
particle production from a color field, formation of a 
QGP and subsequent hadronization. Correlation and 
flow studies suggest that the lifetime of the system is 
short (< 10 fm/c) and, for the first time, there is evi- 
dence suggesting thermodynamic equilibrium already at 
the partonic level." 



V. ANNOUNCEMENT OF PERFECT LIQUID 

BNL and other RHIC-related institutions announced 
in 2005 formation of a "perfect liquid" in high-energy 
Au-Au collisions [115] . The term refers generally to a 
fluid with very low shear viscosity 77, perhaps approaching 
a theoretical quantum lower limit denoted by 77/s = 
(natural units), where s represents entropy density. The 
QCD fluid reportedly formed in Au-Au collisions is de- 
scribed as a strongly- coupled QGP or sQGP, in contrast to 
the weakly-coupled QGP (ideal gas of quarks and gluons) 
anticipated at large energy densities and temperatures 
corresponding to QCD asymptotic freedom. Theoretical 
arguments favoring perfect liquid formation focus on in- 
ference from data of strong hydrodynamic flows and the 
systematic properties of those flows, especially elliptic 
flow represented by quantity i>2. Low viscosity implies 
hydrodynamic expansion with little dissipation. Addi- 
tional evidence for sQGP is sought in modifications to 
QCD jet formation and bulk properties of the medium. 

This section reviews specific theoretical arguments pre- 
sented in Refs. [2Tjm] as representative. In Ref. [5T] it is 
acknowledge that interpretation of particle data can be 
ambiguous: "...a healthy bit of 'luck' was essential in or- 
der to find the 'needles in the haystack' that are least dis- 
torted by uncertain non-equilibrium hadronic final-state 
dynamics." The conjectured 200 GeV collision scenario 
is (i) interpenetration of two thin sheets of "color glass" 
(gluon condensate within projectiles Lorentz contracted 
to 0.1 fm) (ii) initial ID (longitudinal) Hubble expan- 
sion of collision products, (iii) produced quarks and glu- 
ons thermalize locally, (iv) the resulting QGP develops 
collective flows described by hydrodynamics, (v) quarks 



and gluons finally hadronize, (vi) hadrons may continue 
to scatter until kinetic decoupling and free streaming. 

The initial conditions (e.g., gluon and energy densities) 
are inferred from the final-state hadron multiplicity and 
transverse energy E t . Extrapolation backward from the 
final state assumes the system density falls as 1/t (proper 
time). Expansion without work (no secondary scattering 
of particles) would maintain a constant energy per parti- 
cle. In a thcrmalized system entropy is conserved (isen- 
tropic expansion), entropy per particle remains constant 
and the medium cools. Onset of 3D expansion occurs 
shortly before freezeout. That argument combined with 
lattice QCD (LQCD) results provides a basis to claim 
sufficiently high initial energy densities for QGP forma- 
tion. The principal empirical evidence for conjectured 
sQGP follows. 

Collective flow - A claim of thermalized bulk mat- 
ter requires observation of collective flow, the primary 
observable. Flow tests two conditions: (a) thermaliza- 
tion and (b) validity of an equation of state or EoS (a 
relation between energy density and pressure). The ob- 
served flow is said to be consistent with a QCD EoS. 
Elliptic flow V2 is emphasized because it is thought to 
develop early in the collision and to be more sensitive 
to QGP formation. Measured values for central Au-Au 
collisions are compatible with nonviscous ideal hydrody- 
namics below p t fa 2 GeV/c. A QCD EoS is said to be 
confirmed by the extent of mass ordering of V2(pt) for 
identified hadrons below 2 GeV/c. The validity of ideal 
hydro for observed V2 implies very strong interactions of 
quarks and gluons early in the collision and therefore lit- 
tle dissipation. "The smallness of dissipative corrections 
[required for hydro descriptions of «2 data]... is in itself a 
remarkable and unexpected discovery. [...] ...the QGP at 
RHIC is almost a perfect liquid. [...] Elliptic flow mea- 
surements confirm. . .local thermal equilibrium. . ." early in 
the collisions [H]. 

Jet quenching or suppression - A hydrodynamic 
description of V2 data is said to break down at higher p t 
values. For instance, the V2(pt) mass ordering reverses 
above 2 GeV/c. Local equilibrium is not achieved be- 
cause of weaker coupling (and therefore more dissipation) 
at higher momentum scales (QCD asymptotic freedom). 
In the interval 2-5 GeV/c the applicable theory is uncer- 
tain. Above that interval perturbative QCD (pQCD) is 
applied to describe jet production and possible jet mod- 
ifications (quenching) within the sQGP. 

The second major discovery at RHIC is the strong sup- 
pression of hadron spectra above 4 GeV/c in more-central 
Au-Au collisions. Since jet fragments should dominate 
that pt interval the suppression effect is referred to as "jet 
quenching." The principal mechanism is thought to be 
parton energy loss via gluon bremsstrahlung within the 
dense colored medium. Jet quenching in turn can be used 
to study the gluon density profile dN g /dy of bulk QCD 
matter. The value inferred from jet quenching is said to 
be compatible with that inferred from (a) the final-state 
hadron density, (b) the initial conditions required for hy- 
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dro descriptions of elliptic flow and (c) estimated CGC 
(gluon saturation) initial conditions based on a saturation 
scale Q s value 1-1.5 GeV. The CGC inferred gluon den- 
sity is much greater than that predicted by the minijet- 
based HIJING Monte Carlo with a default parton spec- 
trum lower bound po = 2 GeV (parton energy) that lim- 
its minijet production to less than 1000. One should also 
note the experimental lower bound Qo ~ 6 GeV (dijet 
energy) for dijet production in p-p collisions [4"S1 146) . 

Jet-related angular correlations may provide tomo- 
graphic information about the dense medium. In p-p 
and peripheral Au-Au collisions in-vacuum dijet ID az- 
imuth correlations include two peaks at and n. In cen- 
tral Au-Au collisions the so-called "away-side" peak at 
7r is greatly reduced or disappears (for certain p t cuts). 
The result is termed a "monojet" and is interpreted to 
demonstrate that one jet of a triggered dijet is essen- 
tially absorbed by the dense medium. The away-side 
peak is restored in d-Au collisions, demonstrating that 
cold nuclear matter in a gold nucleus does not reduce 
the initial-state gluon density by "shadowing." "The 
observed jet quenching in Au-Au [collisions] is due to 
parton energy loss." "Theoretical analysis of jet quench- 
ing... strengthens the case for multiple strong interactions 
of the quark and gluon constituents of the matter made 
at RHIC." 

CGC initial conditions - The collision initial con- 
ditions (IC) must be specified for any hydro description 
of flow data. The validity of hydro descriptions and in- 
ference of QGP formation is said to depend on the accu- 
racy of IC estimates as well as a theoretical EoS. Mini- 
jet production (as modeled by HIJING) and the CGC 
prediction represent alternative IC descriptions. Initial 
comparisons of predictions for hadron production by the 
competing models with data led to preference for the 
CGC description, "...the surprising very weak centrality 
and beam energy dependence [slow growth] observed [in 
the data] is most satisfactorily explained and predicted 
by the CGC..." The reference is to HIJING predictions of 
more-rapid increases. The comparison, seen to rule out 
minijets and the TCM in favor of the CGC, "...is one of 
the strongest lines of empirical evidence..." for the CGC 
IC (see comments in Sec. XI B I. 

Conclusions - Criteria for QGP discovery: (a) Nu- 
clear matter is created at large enough energy densities 
that it must consist of quarks and gluons when compared 
with LQCD. (b) The matter must be thermalized. (c) 
Properties of the bulk matter must be consistent with 
QCD predictions based on hydro, LQCD and pQCD. 
The published RHIC data are said to satisfy those re- 
quirements. Ideal hydro agrees to a surprising degree 
with some V2 data, implying that the matter must be a 
strongly-coupled QGP. Other evidence from data seems 
to indicate a prominent role for "constituent" or "va- 
lence" quarks in flow manifestations [constituent-quark 
scaling of Vsipt) data] and spectrum structure (evidence 
for quark recombination in PID spectrum systematics). 

Theoretical arguments for a strongly-coupled QGP or 



"perfect liquid" formed in RHIC Au-Au collisions were 
based on data from the initial RHIC running periods. 
Much work has been carried out since to test theoretical 
conjectures and provide a more-detailed view of collision 
dynamics and the nuclear matter created at RHIC. We 
review more-recent results for a subset of topics that seem 
most closely related to sQGP claims. 



VI. HADRON YIELDS AND SPECTRA 

Particle production mechanisms, bulk matter thermo- 
dynamics at low pt including radial and elliptic flow, 
unexpected hadron ratios at intermediate pt and jet 
quenching at high p t continue to be explored. With in- 
creasing beam luminosity and upgraded detector compo- 
nents identified-particle (PID) spectra for more-massive 
hadrons and extension to higher p t greatly increase the 
available spectrum information in the three p t intervals. 



A. Hadron yields near mid-rapidity 

As noted, the systematics of hadron yields within some 
fixed angular acceptance (A?7, A(f>) and the hadron den- 
sity distribution on r\ are essential to establish the IC for 
hydro calculations and to test conjectures about hadron 
production mechanisms (e.g., soft vs hard particle pro- 
duction in a two-component model). 

Figure 11 (left) shows a comparison of participant- 
scaled hadron yields within |?7| < 1 vs Au-Au central- 
ity (points) for 20 and 200 GeV collisions 47 . The 
data are compared with three theoretical models: (a) 
HIJING minijet two-component model (dashed curve), 
(b) fitted two-component model (dotted curve) and (c) 
saturation-scale (CGC) model (solid curve). It is impor- 
tant to note that the data span only the most-central 
40% of the Au-Au cross section. It is concluded that the 
CGC model best describes the data: 

"We find that the Hijing calculation gives the expected 
increase of pQCD minijet production with centrality over 
this energy range, but the predicted increase is now in 
strong contradiction to the data. The flat centrality de- 
pendence of the [energy] ratio is relatively well described 
by the parton saturation [CGC] model calculation" [47] . 

Figure 11 (right) shows the same data normalized by 
the p-p values. The trends for two energies are very sim- 
ilar. The inset compares the HIJING prediction to the 
same ratios. There appears to be a dramatic difference 
that falsifies HIJING. If the data in the left panel are 
taken in ratio 200/20 the result is consistent with a con- 
stant value independent of centrality: "The ratio of the 
measured yields at 200 and 19.6 GeV shows a clear geom- 
etry scaling over the central 40% inelastic cross section 
and averages to R200/19.6 = 2.030.02(stat)0.05(syst). A 
large increase in yield from hard processes, which con- 
tribute to multiplicity, is not apparent in the data, even 
over an order of magnitude range of collision energy." 
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FIG. 11: Left: Charge multiplicity density per participant 
pair vs N pa rt for 20 and 200 GeV Au-Au collisions (points) 
compared to three theoretical predictions (curves) [47]. Right: 
Data in the left panel scaled by p-p charge density showing 
similar trends with centrality for the two energies. 



In Ref. [48] Figs. 13 and 14 extend the above study to 
Cu-Cu, d-Au and several energies. Relative and absolute 
variations of soft and hard components with collision en- 



ergy are discussed in Sec. XI F of the present paper 



B. Low-pt spectra and bulk medium properties 

Hadrons from the lowest pt interval 0-2 GeV/c are 
expected to emerge from a flowing bulk medium. The 
corresponding spectrum shape systematics should re- 
veal medium properties and flow parameters. Identified- 
hadron abundances should test statistical-model predic- 
tions and medium properties such as chemical and kinetic 
freezeout temperatures T c ^ em and T^ in . Certain "scaling 
relations" are sought as an indicator of underlying sim- 
plicity in the sQGP. 

Figure 12 (upper) shows per-participant yield ratios for 
three hadron species, comparing Au-Au and Cu-Cu colli- 
sions on participant number N part |49j . The parameters 
are obtained from blast-wave (BW) fits to PID spectra 
for pions, kaons and protons over a limited low-p t inter- 
val < 1 GeV/c. The yield data for the two collision sys- 
tems do not coincide when plotted on participant number 
Np ar t, do not scale. 



Figure 12 (lower) shows apparent scaling when BW 
parameters are plotted on total hadron multiplicity 
dn c h/dr]. The first panel is T^m vs (/?t), the second shows 
Tchem (upper) and Tki n (lower) vs multiplicity, the third 
shows (/3t) vs multiplicity. The trends for Cu-Cu and 
Au-Au appear to coincide within overlapping multiplic- 
ity intervals. Those results are interpreted to indicate 
that a thermalized bulk-medium description is valid, and 
that "...the kinetic freeze-out properties are determined 
by the initial state." 

Reference [IS] concludes: "These multidimensional 
systematic studies reveal remarkable similarities between 
the different colliding systems. [...] The bulk properties 
studied have a strong correspondence with the total par- 
ticle yield. [...] Within thermal models this reflects a re- 
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FIG. 12: Upper: Per-participant identified-hadron produc- 
tion relative to p-p vs A-A centrality measured by N pa rt for 
Au-Au (light points) and Cu-Cu (dark points) collisions [49] . 
Lower: Bulk- medium properties Tfe in , T c h em and (Pt) vs cen- 
trality measured by dn c h/dr] for 200 GeV Au-Au and Cu-Cu 
collisions. 



lation between the energy per particle at freeze-out and 
the entropy derived from particle yields, which reflects 
the initial state properties for adiabatic expansion. [...] 
The lack of scaling for N part "...suggests that Npart does 
not reflect the initial state of the system accurately." 



C. Intermediate pt and the baryon/meson anomaly 

Hadrons from the intermediate-pt region 2-5 GeV/c 
may arise from thermal hadron production, parton frag- 
mentation to jets or a hybrid combination. The possi- 
bility that hadrons are formed by "constituent quark" 
recombination or coalescence from a thermalized QGP 
has received much theoretical attention [5"UH5"2"] . The 
so-called "baryon-meson" puzzle is intriguing. In more- 
central Au-Au collisions the baryon-to-meson ratio de- 
parts substantially from expectations for fragmentation 
to jets as in elementary collisions, increasing from « 0.3 
near 10 GeV/c to exceed unity at its maximum value 
near 2.5 GeV/c. 

Figure [13] shows the proton-pion ratio plotted vs 
hadron p t for particles (left) and antiparticles (right) 
and for central (upper, black points) peripheral (lower, 
red points) and d-Au collisions (open points) [53] . The 
dash-dotted curve in the left panel represents a prediction 
from fragmentation functions [54) . The dashed and dot- 
ted curves represent predictions from two recombination 
theories [50] [51] . Within the intermediate p t region the 
theories seem to describe the central Au-Au data well, 
implying that the ratio anomaly in central collisions may 
signal hadron formation from a QGP in which constituent 
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pared to a fragmentation function parametrization derived 
from e + -e~ collisions (dash-dotted curve) and to recombina- 
tion theories (dashed and dotted curves) . 



quarks play a significant role. 



D. High-pt spectra and jet suppression 

Hadrons from the high-p t region above 5 GeV/c should 
be dominated by parton fragmentation to jets. The 
principal issue is possible jet quenching in more-central 
A-A collisions by parton energy loss in the dense col- 
ored medium or sQGP. How does jet quenching change 
with system size, both A-A centrality and A? How does 
jet quenching depend on the leading parton (quark or 
gluon)? 



Figure 14 (left) shows Raa measurements for central 
and peripheral Au-Au and several centralities of Cu-Cu 
collisions extending to 8 GeV/c [55 . The suppression 
in most-central Cu-Cu is substantially less than that for 
Au-Au, which is not surprising for the smaller collision 
system. Figure 14 (right) shows Raa mean values within 
5-8 GeV/c plotted on N part to reveal a common centrality 
trend for the two collision systems. A similar pair of plots 
(not shown here) demonstrates that the mean p/ir ratio 
within the interval 3-4 GeV/c shows a similar universality 
when plotted on N part . 

Reference [SS] concludes: "The [Cu-Cu] data are found 
to exhibit similar systematic trends over a wide range of 
transverse momenta as Au-Au collisions at the same en- 
ergy with a similar number of participants. [...] The 
participant coverage in these Cu-Cu collisions is in a re- 
gion where the [parton] suppression effects are turning 
on. [...] A detailed study of the proton to pion ratio 
reveals similar systematic dependencies to that found in 
Au-Au data. [...] ...these results indicate similar partonic 
energy loss for both gluons and quarks. [...] ...the [par- 
ton] suppression for different collision species is found to 
be invariant for the same number of participants." 
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FIG. 14: Left: Spectrum ratio -Raa(pj) for several cen- 
tralities from 200 GeV Au-Au (blue) and Cu-Cu (red) col- 
lisions 55.. Right: Raa averaged over pt interval 5-8 GeV/c 
vs Np ar t demonstrating apparent participant scaling of jet 
suppression. 
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| introduces an alternative suppression 

in which the factor l/N^n is replaced 
in the spectrum ratio by l/N part . Whereas the conven- 
tional ratio Raa shows strongly increasing suppression 
with increasing centrality and energy above 50% frac- 
tional cross section the Raa^ trend is nearly indepen- 
dent of Au-Au centrality. Reference [56] remarks: "The 
apparent dominance of the initial-state geometry [-/V par . t ], 
even for observables closely related to the dynamical evo- 
lution of heavy-ion collisions, is one of the key features 
of these interactions that remains to be understood." 

Reference (57) compares PID particle spectra and ra- 
tios in p-p and Au-Au collisions, extending measurements 
up to 15 GeV/c. Variation of particle-to-antiparticle ra- 
tios confirms the expected changing mixture of quark 
vs gluon jets with increasing p t (as the x of the hard- 
scattered parton increases toward valence quarks). The 
available theoretical parametrizations of PID fragmenta- 
tion functions describe the pion data but are not con- 
sistent with the kaon and proton data. Several hadron 
species share a common Raa suppression above 6 GeV/c. 

Some studies have employed theoretical models of par- 
ton energy-loss in the sQGP to estimate the color-charge 
density. Reference [58] uses Raa angular variations rel- 
ative to the reaction plane to infer the pathlength de- 
pendence of jet suppression. Greater suppression is ob- 
served out of plane than in plane. Suppression is ap- 
proximately independent of pt above 4 GeV/c and varies 

2/3 

with centrality approximately as N p ' art consistent with 
theoretical models. The inferred path-length dependence 
appears to conflict with pQCD energy loss theory. Ref- 
erence [59] describes an ambitious program to establish 
a self-consistent description of Raa data with global fits 
by energy-loss models. By relating the fit parameters to 
pQCD energy-loss theory the goal is to infer the system- 
atics of the color-charge density in the sQGP. 

These high-p t studies emphasize jet suppression as in- 
ferred from spectrum ratio Raa over a pt interval above 
4 GeV/c. That interval includes only a small minority of 
the hadron fragments within any jet [41, 145) . From such 
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analysis we do not learn what happens to the energy ap- 
parently "lost" from the leading parton, or whether it is 
lost at all from an intact but modified jet. Methods that 
access jet structure at lower p t are required to pursue 
such questions. 



VII. v 2 AND ELLIPTIC FLOW 

Because of its importance in providing the primary ev- 
idence for a strongly- coupled QGP with apparently very 
low viscosity elliptic flow vi has continued to receive 
a great deal of attention, with a growing number of 
subtopics. 



A. V2 centrality, n and energy dependence 

The initial 130 GeV Au-Au p t -integral V2 data leading 
to inference of "perfect liquid" were extended to 20, 62 
and 200 GeV, to the smaller Cu-Cu collision system and 
to a variety of hadron species. In Ref. [60] the strongly- 
peaked r\ dependence of i>2 over a large r\ acceptance was 
confirmed. The T) dependence does not change signifi- 
cantly with centrality within the upper 50% of the Au-Au 
total cross section. The strong n dependence is a signifi- 
cant problem for hydro models of elliptic flow. Collision 
energy dependence is also of interest. Does u 2 continue 
its rise from Bevalac-AGS to 200 GeV? Does it saturate 
at higher energies? 



B. Nonflow and flow fluctuations 

v-i data may include a systematic error or bias from 
"nonflow" defined as contributions to v-i not related to 
the A-A reaction plane |23j . A number of strategies or 
"methods" denoted by V2{method} have been developed 
in attempts to reduce or eliminate nonflow. w 2 data may 
also include contributions from non-Poisson fluctuations 
possibly corresponding to initial-state A-A geometry fluc- 
tuations. The two issues are intimately related. 

The origin of nonflow remains ambiguous because by 
hypothesis jets should not contribute to hadron produc- 
tion below 2 GeV/c where flows are thought to dominate, 
are therefore largely discounted as a nonflow mechanism. 
More-recent V2 methods attempt to reduce nonflow bias. 

Figure [l5| (left) shows i>2(6) data for several analysis 
methods |23j . "Standard" denotes the original event- 
plane (EP) method [ST]. ^{2} and t>2{4} respectively 
denote two- and four-particle cumulant methods. Other 
methods represent specific pair acceptances on (?7i , 772) 
and other numerical algorithms (e.g., scalar product, Q- 
vector). The data from different 1*2 methods show sub- 
stantial differences, mainly for most-central and most- 
peripheral collisions, that may be due to nonflow or flow 
fluctuations. Such differences are used to infer system- 



atic uncertainties due to nonflow bias, but the various 
methods share certain assumptions in common. 
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FIG. 15: Left: pt-integral v% vs fractional cross section from 
200 GeV Au-Au collisions for several V2 methods [23 . The 
differences among methods are conventionally invoked to es- 
timate systematic uncertainties due to nonflow. Right: Com- 
parison of «2 and 174 for 200 GeV Au-Au collisions |62) . The 
close agreement among the employed methods suggests min- 
imal nonflow bias for those methods. 



Figure 15 (right) shows V2 (elliptic) and V4 (hexade- 
capole) flow data vs N part for several pt bins 62 . The 
various plotting symbols represent several V2 methods 
in which different r\ acceptances (detector elements) are 
used to estimate the event plane. The close agreement of 
the different methods is interpreted to imply that non- 
flow bias is made negligible by EP estimates from large-?7 
detectors [53] • It is notable that the large-77 results agree 
closely with w 2 {2} data. 

Elliptic flow is believed to be a final-state manifesta- 
tion of the initial-state A-A overlap region in non-central 
collisions. V2 data are therefore compared with estimates 
of the initial-state geometry in the form of eccentric- 
ity e. Earlier estimates of e were based on an optical 
Glauber model. Reference [64] introduces the concept 
of a Monte Carlo Glauber or participant eccentricity to 
model the initial-state A-A geometry. Whereas the op- 
tical (or standard) Glauber eccentricity falls toward zero 
for central and peripheral A-A collisions the MC eccen- 
tricity remains significantly nonzero for central collisions 
and rises asymptotically toward unity for peripheral col- 
lisions, a dramatically different scenario. The resulting 
difference in i^/e trends has major implications for hy- 
drodynamic interpretations. 

Figure 16 (left) shows a comparison between optical 
(standard)? and Monte Carlo (participant) e vs N part 
for Au-Au and Cu-Cu collisions 64 . In that paper it 
is observed that i^/e for both Au-Au and Cu-Cu follow 
the same trend vs centrality (N part ) if MC e is used, sug- 
gesting that elliptic flow in Cu-Cu is comparable to that 
in Au-Au for the same number of participant nucleons 
even though the Cu-Cu system is much smaller. The 
importance of participant eccentricity (and a participant 
event-plane estimate) in such comparisons is emphasized. 

Based on method definitions and certain assumptions 
the difference w 2 {2} — ^ll^} is said to have two main 
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FIG. 16: Left: A-A overlap eccentricity e vs N part for optical 
(standard, solid points) and Monte Carlo (participant, open 
points) Glauber model simulations and for Au-Au and Cu-Cu 
collisions [64] . Note the offset zero. Right: Estimated final- 
state V2 fluctuations measured by R v (solid points) compared 
to theory estimates of initial-state eccentricity fluctuations 
measured by i? e (open points) vs Au-Au centrality measured 
by dn c h/dri [65] . 



contributions: (a) a nonflow contribution 5 2 from sev- 
eral possible sources and (b) a v% fluctuation contribu- 
tion denoted by 2<r 2 2 [55]. Various attempts are made 
to distinguish the two contributions and to estimate the 
r.m.s. measure of relative v-i fluctuations R v = <7 V2 /v2 in 
comparison with theoretical predictions of IS geometry 
fluctuations in the form R € = er e /e. 

16 (right) shows comparisons between in- 
fluctuations (solid points) and theory (open 



Figure 
ferred V2 
points) [65 



The experimental data represent upper lim- 
its assuming nonflow <5 2 is zero in all cases. However, 
there is nothing to rule out the other extreme, that a% 2 
may be negligible in all cases. Alternative information on 
possible nonflow contributions to V2 (e.g., ^2 arising from 
jets) is available from other angular correlation measure- 
ments I66j . In Ref. |67j the nonflow contribution was 
estimated with (771, 772) correlations. It was assumed that 
the (long-range) nonflow contribution is approximately 
the same as for p-p collisions. That is a strong assump- 
tion since jet contributions may scale as Nun ~ 1000 in 
central Au-Au relative to p-p collisions. 



C. PID V2(pt) and constituent-quark scaling 

Analysis of V2(pt, b) for identified hadrons (particle ID 
or PID) confronts two main issues: so-called mass scaling 
below 2 GeV/c interpreted to confirm conjectured hy- 
drodynamic flow and so-called constituent-quark scaling 
within 2-5 GeV/c interpreted to indicate that the flowing 
medium consists of quarks and gluons. Above 2 GeV/c 
t>2 is said to be driven more by quark content than by 
hadron mass, due to hadron formation by quark coales- 
cence [51] , Quark- number scaling may indicate that in 
RHIC heavy ion collisions collective motion is established 
among quarks and gluons before hadrons are formed, ap- 
parently confirming that a strongly-coupled QGP with 



partonic degrees of freedom is created in RHIC collisions. 
Experimental questions include the accuracy or signifi- 
cance of inferred scaling trends and the dependence on 
centrality, system size and collision energy. 

Reference [68] reports comprehensive scaling with col- 
lision geometry, system size and transverse kinetic en- 
ergy KE t = m t — rrih- It is reasoned that early pressure 
gradients drive KE t which should then replace pt as a 
plotting variable. At smaller p t or KE t mass scaling is 
observed, leading to inference of hydrodynamic flow. At 
larger KE t mesons and baryons scale separately. If quark 
number n q scaling (of w 2 and KE t ) is included all hadron 
species seem to scale together over a large pt interval. 
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FIG. 17: Left: V2 for identified hadrons vs pt from minimum- 
bias 200 GeV Au-Au collisions, both quantities scaled by 
constituent-quark number n q — 2 or 3 |68) . Right: The same 
data plotted vs transverse kinetic energy KE t = m t —mh also 
scaled by n q showing closer scaling. 



Figure 17 shows constituent-quark scaling on p t (left) 
and KE t (nght) for minimum-bias 200 GeV Au-Au data. 
In the right panel several hadron species appear to follow 
the same trend (true for ir, K, p, A, 5). Scaling universal- 
ity is also reported between Cu-Cu and Au-Au systems 
and for various centralities of each. Those data are ob- 
served to scale (follow a universal functional form) when 
plotted as V2(pt,b)/v2(b) vs pt- A more recent study re- 
ports that strong deviations from such scaling in noncen- 
tral Au-Au collisions (less than 20% central) are observed 
at larger KE t values (KE t /n q > 1.5 MeV/c 2 ) [53] • 

Reference [7D] presents a comprehensive study of PID 
v^iPtib) for several hadron species for p t up to 6 GeV/c 
and a range of Au-Au centralities. Figure 18 shows PID 
V2 data for 200 GeV Au-Au collisions and several cen- 
trality intervals [70]. The figure illustrates agreement 
between ideal hydro (curves) and data (points) below 2 
GeV/c in mass ordering and magnitude, seeming to sup- 
port a hydrodynamic interpretation. It also shows sepa- 
ration of meson and baryon data in the intermediate pt 
interval 2-5 GeV/c interpreted to indicate constituent- 
quark scaling. 

In Figs. 9-11 of Ref. [70] a differential study of devia- 
tions from a common scaling trend are shown. Substan- 
tial systematic deviations are observed, both in data and 
in theory predictions, leading to the conclusion that nei- 
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FIG. 18: V2 vs m t — mi, for identified hadrons and several 
ccntrality intervals from 200 GeV Au-Au collisions (points) 
compared with ideal hydrodynamic predictions (curves) [TO] . 



ing within each centrality bin but not between centrality 
bins. The curves at upper right are hydro predictions 
that approximately follow the data (with mass scaling) 
below 2 GeV/c but strongly overshoot above that point. 
Hydro theory shows no quark-number scaling. It is con- 
cluded that "... mass ordering at low p t alone is not 
sufficient to claim thermalization in Au-Au collisions at 
RHIC." 

When the same data are also scaled by p t -integral (^2) 
for each hadron species the scaling at larger p t seems to 
improve as in panel (c) . When the scaling is done instead 
with {v2)ch common to all charged hadrons (mainly pi- 
ons) the scaling improves further as in panel (d). In 
the same panel (inset) significantly negative V2 values 
are shown for the first time, interpreted to support the 
strong bulk expansion inferred from spectrum analysis, 
but see Fig. 35 (right). 



ther V2 data nor hydro theory follow strict mass scaling 
or constituent-quark scaling as they are usually defined. 



D. V2 from the RHIC beam energy scan 
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FIG. 19: v^ipt) data plotted vs pt and mt — for several 
hadron species and for several centrality intervals of 200 GeV 
Au-Au collisions scaled by constituent quark number n q and 
by pt-integral (112) for individual species and all hadrons [70] . 
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1 [G. 20: Upper: V2{i}(p t ) vs p t for several collision energies 
showing similar trends for different energies [71] . Significant 
variation apparent at lower-pt is compared to no variation at 
larger p t . Lower: 112 vs r\ for several collision energies showing 
the same functional form on r)/ybeam for all energies. 



Figure [19] shows constituent-quark scaling systematics 
for K, A and S hadrons [70]. In panel (a) there is no scal- 
ing with participant eccentricity for p t > 2 GeV/c. For 
p t < 2 GeV/c there is mass scaling with m t — rrih, and in 
the intermediate region there is constituent-quark scal- 



The RHIC beam energy scan is intended to determine 
the energy where sQGP production "turns off," and pos- 
sibly locate a conjectured "critical point" on the QCD 
phase diagram. Figure 20 (upper) shows first v 2 results 



for 7.7-39 GeV from STAR [71]. The BES data are com- 
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pared with other V2(j>t) data up to Pb-Pb collisions at 
2.76 TeV. There is surprisingly little variation: small 
monotonic increase at the lowest p t but no significant 
increase in the intermediate-p t region. 



Figure 20 (lower) shows evolution of the rj dependence 
of V2 over the range 7.7 - 200 GeV indicating substan- 
tial variation. However, if the ^2(^7) trend is rescaled by 
the beam rapidity the functional form appears to be uni- 
versal (right panels). It is also noted that the difference 
w 2{2} ~ conventionally attributed to a combina- 

tion of nonflow and flow fluctuations is reduced at lower 
energies. Viscous hydrodynamics does not reproduce the 
energy dependence of v^ipt). 



E. Higher harmonic flows 



Measurement of higher harmonic flows (t>3, W4, etc.) 
in heavy ion collisions is an outgrowth of attempts to 
measure flow fluctuations as a response to initial-state 
geometry fluctuations [55] and replacement of the con- 
jectured Mach-cone interpretation of away-side ID peak 
distortions (as inferred from ZYAM subtraction) with 
conjectured "triangular flow" represented by V3 as one 
of several higher harmonic flows [72] [73] . 

Reference |74j reports measurements of 1)2, v 3 , V4 vs pt 
and N part for 200 GeV Au-Au collisions. The data are 
apparently well-described by hydro models incorporating 
a Glauber model of initial-state geometry with fluctua- 
tions below 2 GeV/c. The agreement appears to provide 
evidence for such geometry fluctuations. The results are 
interpreted to support a value for shear-viscosity mea- 
sure 77/s rs 1/47T, the quantum lower limit consistent with 
"perfect liquid." 

Reference [75] reports extraction of u|{2}(A?7) by 
Fourier decomposition on azimut h o f slices bin-by-bin 
across the eta acceptance. Figure [2~i"| (left) shows vf{2} 
extracted for three Au-Au centrality bins and for CI, LS 
and US charge combinations. Narrow and broad peak 
structures are apparent. Those results can be compared 
directly with 2D angular correlations presented in Fig. 31 
(right) from Ref. [5j. 



Figure 21 (right) shows centrality trends for several V3 
methods. For some methods there is substantial trian- 
gular flow inferred even for very peripheral Au-Au colli- 
sions. It is notable that the v 3 values inferred from data 
in the left panel coincide with those labeled v%{TPC} 
(event-plane determined within the TPC 77 acceptance, 
right topmost data) that describe all angular correlations 
including jets. A substantial role for nonflow is rejected 
because LS and US results are similar. Based on model 
comparisons it is concluded that the W3 structure "...is 
mainly due to Ar) dependent fluctuations" [75] . 
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FIG. 21: Left: Triangular flow measured by fl{2} inferred 
from Fourier fits to ID projections of 2D angular correlations 
onto ID azimuth from narrow Arj bins [75] ■ Right: Compar- 
ison of v$ inferred from various methods. Data from the left 
panel (wide Gaussian) correspond to t>3{2} inferred within 
the STAR TPC acceptance (upper-most points). 



VIII. JET ANGULAR CORRELATIONS 

Several methods are employed to infer jet angular 
structure in A-A collisions, including triggered ID az- 
imuth correlations with background (ZYAM) subtrac- 
tion, triggered 2D angular correlations on (77, (f>), event- 
wise jet reconstruction and untriggered (minimum-bias) 
combinatoric ID and 2D angular correlations. 



A. ID azimuth correlations and ZYAM subtraction 

Triggered dihadron correlations on azimuth are formed 
by selecting the highest-p t trigger particle in each event 
and forming the pair distribution on azimuth difference 
A<p = <t> — 4>trig relative to that particle. A subset of the 
remaining particles in the event is referred to as asso- 
ciated. Conditions are placed on accepted trigger and 
associated p t ranges. The trigger particle is assumed 
to act as a proxy for the parent parton of a (triggered) 
jet. The result is jet angular correlations plus a com- 
binatoric background including a sinusoid contribution 
interpreted as elliptic flow. The background is estimated 
using measured vi data relevant to collision and trigger 
conditions. The overall background amplitude is deter- 
mined by the zero-yield-at-minimum (ZYAM) condition, 
which assumes that jet-related peaks do not overlap |76j . 

Figure [22] compares ZYAM-subtracted trigger- 
associated angular correlations for trigger particles 
in 4-6 GeV/c from 200 GeV central Au-Au and p-p 
collisions [77]. The upper panels include associated 
particles within 0.15-4 GeV/c, the lower panels include 
only 2-4 GeV/c. In the former all jet-related structure 
increases from p-p to central Au-Au and the AS peak is 
undistorted. In central Au-Au collisions it is observed 
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FIG. 22: Jet-related dihadron correlations on A-q and A(f> 
inferred from 0-5% central 200 GeV Au-Au and p-p collisions 
by ZYAM subtraction and trigger-associated pt cuts [77] . 



associated p t cut combinations |H0]- For the lower trig- 
ger cut 3-4 GeV/c AS double peaks are observed near 
A<fi = 7r ± 1.1. The double peaks are said to appear in 
a "shoulder" region (SR). Jet structure in more-central 
Au-Au is interpreted to have four components: (i) a jet 
fragmentation component near Acj> = 0, (ii) a punch- 
through jet fragmentation component near A(f> = 7r, (iii) 
a medium-induced component near A<p = and (iv) a 
medium-induced components near A<p = 7r ± 1.1. The 
AS double-peak separation is said to be independent of 
p t as expected for conjectured Mach shocks (Mach cones). 
Similar structures were reported by Ref. [81] where AS 
peak evolution was described as follows: "The transition 
from a broad away-side structure at low pt to a narrow 
structure at higher p t would then signal the change from 
away-side structures dominated by bulk particle produc- 
tion from the medium to a situation where jet-fragments 
dominate." 



that the away-side peak includes a softer particle com- 
plement than for p-p collisions. It is noted that the AS 
peak in the upper-left panel has a dipole cos(A</>) shape, 
suggesting global momentum conservation as the main 
source. However, the minimum-bias AS peak for p-p 
collisions has a similar shape [30] EH [78] . It is concluded 
that jets are modified in more-central Au-Au collisions. 

Further exploration of dihadron correlation systemat- 
ics with various trigger-associated cut combinations ap- 
peared to reveal strong distortion of the AS peak shape 
for some combinations. A possible source might be pro- 
duction of Mach cones as a medium response to interac- 
tion of energetic partons with the QCD medium [79] . 
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FIG. 23: Jet-related dihadron azimuth correlations from 
200 GeV Au-Au collisions for various combinations of trigger- 
associated p t cuts [80] . 

Figure [23] compares dihadron correlations from 200 
GeV Au-Au and p-p collisions for several trigger- 



B. 2D angular correlations and the "ridge" 

Trigger particles can also be used to generate 2D an- 
gular correlations on (77, <f>) providing more information 
about jet structure. The systematic variation of 7/ elon- 
gation first observed for minimum-bias correlations re- 
ported in Ref. [B] is then revealed with trigger-associated 
cuts, and a same-side "ridge" structure is identified. 



d+Au minimum bias Au+Au central 




FIG. 24: 2D angular correlations relative to a 3-4 GeV/c 
trigger-particle (angle relative to trigger) for 200 GeV d-Au 
collisions (left) and central Au-Au collisions (right) |82| . 

Figure [24] shows triggered 2D angular correlations from 
d-Au (left) and central Au-Au (right) collisions [52] ■ The 
trigger p t interval is 3-4 GeV/c and associated interval 
is 2-3 GeV/c. In the d-Au case the expected nearly- 
symmetric SS 2D peak attributed to jets is observed. In 
the Au-Au case "long-range" (on rf) structure emerges 
that extends beyond the STAR TPC acceptance. The 
elongated structure is referred to as a "ridge." System- 
atic studies are interpreted to conclude that "jet-like" 
and "ridge-like" structures come from distinct mecha- 
nisms. The latter (presumed nonjet) mechanism has re- 
ceived much theoretical attention [83]. In Ref. [84] it is 
shown that the same-side ridge appears to extend to large 
IA77I « 4. 
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Triggered dihadron correlations may be used to probe 
the QCD medium (jet tomography). Strong jet quench- 
ing and apparent large reduction of AS jet yields suggests 
formation of an opaque colored medium in more-central 
collisions. The parton pathlcngth through the medium 
may then strongly influence jet modification. 

In Ref. [55] direct photons are used to establish a trig- 
ger direction relative to which away-side jet structure 
is studied. The photon energy then estimates the re- 
coil energy of the parent parton of the away-side (part- 
ner) jet. Reduction of the away-side associated-hadron 
yield shows no significant dependence on photon (par- 
ton) energy. It is concluded that dependence on parton 
species (light quark or gluon) and parton pathlength in 
the medium is also small. 

In Ref. [HD a novel rfowWe-trigger plus associated (2+1) 
technique is introduced: two high-p t triggers are sepa- 
rated by 7r on azimuth. The double trigger is said to 
bias toward tangential emission of jet partners outside 
an opaque core. Such dijets then do not interact with 
the core. Double-triggered same-side and away-side jet- 
related peaks in central Au-Au have similar characteris- 
tics to those in d-Au collisions. 



Figure [25] shows a study of the structure of the SS 2D 
peak vs several collision parameters (system size, cen- 
trality and energy) |87j . The upper panel shows SS peak 
shapes on Arj for specific cuts on A(f> illustrating evolu- 
tion of r) elongation with collision conditions. The lower 
panel shows SS peak shapes on A<fi illustrating the A(f> 
cuts and ZYAM subtraction (curves). The SS peak is de- 
scribed as having jet-like (narrow on Arj) and ridge-like 
(broad on Arj) components. The ridge structure is said 
to persist to pt,trig ~ 6 GeV/c and Pt,assoc ~ 3 GeV/c. A 
2D background is subtracted according to the ZYAM pre- 
scription. It is interesting that for 0-12% central Au-Au 
(upper panel, lower right) the elongated "ridge" struc- 
ture is consistent with monotonic decrease toward zero 
on At], possibly part of a single monolithic peak. In gen- 
eral, both jet-like peak and ridge share the same small 
azimuth width. The ridge persists in Cu-Cu collisions 
and at 62 GeV. And the energy dependence of jet-like 
and ridge- like structures is equivalent. 



IX. FLUCTUATIONS 

Fluctuation measurements were expected to reveal ev- 
idence for traversal of the QCD phase boundary in the 
case of a first-order transition. However, fluctuation 
signals are thought to be less important for a smooth 
crossover transition. More recently, the search for a con- 
jectured QCD critical point (CP) in association with the 
RHIC beam energy scan (BES) program has been advo- 
cated. Fluctuations of event- wise mean p t , multiplicity, 
net charge, net baryon number and hadron species ratios 
such as the K/n ratio are reported. Several fluctuation 
analyses are summarized below. 

Reference [55j describes (pt) fluctuations measured by 
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FIG. 25: Systematic study of same-side 2D peak shapes for 
different collision systems [ST]. Upper: Peak shapes on Arj 
for specific A<f) acceptance cuts and Cu-Cu, d-Au and Au-Au 
collisions at 62 and 200 GeV. Lower: Peak shapes on A(j> for 
the same systems showing ZYAM subtraction. 



(Apt^Apt.j), where Ap t — p t — pt is the particle p t de- 
viation from the ensemble mean. The definition leads to 
a per-pair fluctuation measure. The general trend with 
Au-Au centrality is decrease approximately as 1/N part 
or l/n c h as expected, and the increase with energy is 
quite slow. If that measure is multiplied by dn c h/dr] 
(converting to a per-particle measure) the centrality de- 
pendence is similar to Fig. 10 (left) from Ref. [Mj and 
there is a strong increase with collision energy. The pa- 
per concludes that "...there are clear nonzero pt correla- 
tions. [...] The centrality dependence... may show signs 
of...thcrmalization, the onset of jet suppression, the sat- 
uration of transverse expansion in central collisions, or 
other processes." 

Reference [59] addresses charged-particle density cor- 
relations on r) by measuring multiplicity fluctuations. 
An attempt is made to locate the QCD phase bound- 
ary in the context of the Ginzburg-Landau description of 
critical phenomena. Multiplicity distributions in r\ win- 
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dows of increasing size 5rj and various Au-Au centrali- 
ties are modeled with a negative binomial distribution 
(NBD). The NBD model parameter k measures devia- 
tions from Poisson fluctuations (no correlations) and is 
parametrized by 1/k — 2a^/5r] + (3. 1/k trends on 5n 
and N part represented by parameters a£ and /3 tend to 
decrease as 1/N part since 1/k is also a per-pair measure. 
The study concludes "...The behavior [of product a£] 
may be explained by the onset of a mixture of different 
types of particle production mechanisms which are not 
necessarily related to temperature or density correlations. 
However, interpreted within the Ginzburg-Landau frame- 
work the local maximum of the at; product [on N part , at 
most a two-sigma effect] could be an indication of a crit- 
ical phase boundary." 

Reference [SU] measures fluctuations with so-called 
normalized variance a^ lch /n c h denoted by w c h, a per- 
particle measure. The relation to NBD model parameter 
k is w c h = 1 + n c h/k- An estimate w c h,dyn of the excess 
fluctuations relative to a statistical reference is made. 
The trend of w c h.dyn is similar to that in Fig. u0\ (left) 
from Ref. [35] describing (p t ) fluctuations. The study 
concludes "...there is no evidence of critical behavior re- 
lated to the compressibility observable in this dataset. 
There is also no significant evidence of dynamical fluctu- 
ations that are dependent on the transverse momentum 
or the charge of the particles measured. [...] Although 
this analysis does not observe evidence of critical behav- 
ior, it does not rule out the existence of a QCD critical 
point." 

Reference [HH reports an analysis of K/ir fluctuations 
motivated by the possibility to detect critical fluctuations 
in strangeness enhancement near the QCD phase bound- 
ary. Figure 3 of that reference shows K/tt fluctuations as 
measured by per-pair measure fdy n . The trend is domi- 
nated by the expected l/n c h dependence. In Fig. 4 of the 
reference is plotted (dn c h / dr))vdyn (for charge combina- 
tions LS - green, US - blue and CI - red), a per-particle 
measure which increases slowly with Au-Au centrality. 
The study concludes "These results may indicate that, 
due to later stage hadronic rescattering, the decay prod- 
ucts of resonances are less likely to survive in central 
collisions than in peripheral collisions." 



X. ALTERNATIVE METHODS AND RESULTS 



Section |II B| summarizes alternative analysis methods 
and measures that may be applied without a priori as- 
sumptions about A-A collision mechanisms. In this sec- 
tion we review a selection of those results and consider 
some implications for interpretation of RHIC data. 



A. Mid-rapidity hadron yields 

Mid-rapidity hadron yields are said to play a key role 
in constraining hydro calculations and indicating hadron 



production methods. Two alternative scenarios have 
been proposed: (a) The TCM including longitudinal pro- 
jectile dissociation (soft) and transverse parton scatter- 
ing and fragmentation (hard) and (b) a CGC glasma A-A 
initial state that may transition directly to a QGP. Argu- 
ments have been proposed based on HIJING simulations 
that the TCM and minijets are falsified by hadron yield 
systematica, that the TCM predicts much more yield in- 
crease with A-A centrality than is observed [TT]. In fact 
the TCM with the correct in-vacuum dijet cross section 
(2.5 mb) underpredicts hadron yields for more-central 
Au-Au collisions. The Au-Au data from more-central col- 
lisions requires modified parton fragmentation — increase 
of mean fragment multiplicities by up to a factor 3. 

Figure 26 (left) shows the centrality dependence of di- 
jet production within An = 2 for 200 GeV Au-Au colli- 
sions 92 . The GLS curve is based on a dijet total cross 
section of 2.5 mb consistent with p-p spectrum data |45j . 
The solid curve corresponds to a 50% increase above 
v = 3 inferred from Au-Au spectrum data [1]. 




FIG. 26: Left: Dijet number Uj per Au-Au collision within 
acceptance An vs centrality measure v for 200 GeV Au-Au 
collisions |92) . Binary-collision scaling of p-p dijet cross sec- 
tion 2.5 mb is shown by the dashed curve. Possible increase 
of the dijet cross section for more-central collisions inferred 
from spectrum analysis [4] [45] is describe by the solid curve. 
The upper hatched band indicates the location of the sharp 
transition in minijet characteristics, also where the dijet num- 
ber per unit n in Au-Au significantly exceeds unity. Right: 
The TCM prediction for per-participant hadron production vs 
centrality based on minijet angular correlations (solid curve) 
compared to spectrum data (points) [52]. The dashed curve 
shows a (scaled) CGC prediction [93]. The GLS extrapolation 
(lower hatched band) is based on measured p-p trends 46 . 



Figure 26 (right) shows the total hadron yield in- 
ferred from p t -integral 2D angular correlations (solid 
curve) 92 1 . The volume of the SS 2D peak (number of 
jet-correlated pairs) is combined with the dijet number 
in the left panel to infer the mean fragment number per 
jet n c h.j- The hard-component yield uHaa = Tij n c h,j 
is combined with fixed soft component Snn to predict 
the per-participant-pair hadron yield (solid curve) |92j . 
The points are integrated from spectrum data [4]. The 
dash-dotted line is a fitted two-component model with 
fixed parameter x « 0.1 [8 . The lower hatched region 
represents a GLS extrapolation of p-p systematics, with 
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x « 0.03 corresponding to a fixed 2.5 mb pQCD dijet 
cross section. The yield increase above that value for 
more-central collisions results from the combination of 
a 50% increase in the dijet cross section and a three- 
fold increase in the mean jet fragment multiplicity as in- 
ferred from Au-Au spectrum data. It is notable that the 
mean dijet number for central Au-Au collisions is 30-60 
corresponding to a large jet fragment yield - 30% of all 
hadrons are contained within resolved jets [52] . 

The accurate TCM description of hadron production 
data for all Au-Au centralities is also notable. A CGC 
prediction (dashed curve) has the wrong functional form 
cx \og(8v) and does not predict absolute yields, has been 
scaled to pass through the more-central data points [93] . 
Any successful description of hadron production must ac- 
commodate the interval of A-A transparency below v = 3 
corresponding to 50% of the A-A total cross section [5]- 



B. y t spectra from p-p and Au-Au collisions 

An accurate mathematical model of p-p spectrum 
structure provides an essential reference for A-A colli- 
sions and strongly suggests the underlying hadron pro- 
duction processes. The multiplicity dependence of p t or 
yt spectrum shapes leads directly to a two-component 
model or TCM. No a priori physical model is imposed. 



Figure 27 (left) shows hadron spectra on y t for ten 
multiplicity classes within acceptance |»;| < 0.5 [15]. The 
spectra have been normalized by soft multiplicity n s 
defined iteratively in terms of inferred hard-component 
multiplicity nh and the relation n c h = n s + n^. h c h w 
n c h/2 is the uncorrected multiplicity expressed in terms 
of the corrected multiplicity (corrected for p t acceptance 
and tracking efficiency). The soft-component reference 
So{yt) is the limiting spectrum shape for n c h 0. 




FIG. 27: Left: Spectra for ten n c h classes from 200 GeV p-p 
collisions [46] . Right: Spectrum hard components extracted 
from the spectra in the left panel. 



Figure 27 (right) is the hard component, obtained by 
subtracting So from the spectrum data, described by a 
peaked (Gaussian) distribution Ho(y t ) with amplitude 
approximately cx n c ^ (solid curves). The spectrum sys- 
tematics are summarized (within angular acceptance 2n 



and An — 1) by 

dn ch 



Vtdyt 



= n s S (y t ) + n h H (yt), 



(2) 



where n h /n s cx n ch and H pp = n h H (y t ) [46 



Figure 28 (left) shows the non-single-diffractive (NSD) 



average of the p-p spectrum hard components (points) 
from Fig. 27 The solid curve is a pQCD calculation of 



the p-p spectrum hard component based on a pQCD par- 
ton spectrum bounded below at 3 GeV and integrating 
to 2.5 mb and measured fragmentation functions from 
CDF, LEP and HERA gSHMj- The dash-dotted curve 
is the Gaussian approximation from Ref. |46j . 




FIG. 28: Left: Spectrum hard component representing 200 
GeV NSD p-p collisions (points) [45] . The solid curve is a 
pQCD prediction for the corresponding fragment distribution 
derived from measured fragmentation functions and a dijet 
total cross section of 2.5 mb. Right: Spectrum hard compo- 
nent for 0-12% central 200 GeV Au-Au collisions (points) [4]. 
The solid curve is a pQCD prediction based on a simple mod- 
ification of fragmentation functions 45 . The dotted curve is 
the prediction from the TCM extrapolated from p-p collisions. 



Figure 28 (right) shows the spectrum hard component 
from 0-12%central 200 GeV Au-Au collisions (points) [4]. 
The dashed curve is the p-p Gaussian (with added power- 
law tail) from the left panel. The dotted curve is a GLS 
prediction for central Au-Au (A-A transparency). The 
solid curve is a pQCD description of the central Au-Au 
data based on a single modification of measured frag- 
mentation functions (single-parameter change in a gluon 
splitting function) [45 . Fragment reduction at larger y t 
is balanced by much larger fragment increase at smaller 
y t that conserves the parton energy within resolved jets. 

Figure [29] shows hard-component ratios rAA — 
Haa/H pp for identified pions and protons from five cen- 
trality classes of 200 GeV Au-Au collisions [45] . Whereas 
conventional spectrum ratio Raa includes the spectrum 
soft component that strongly biases the ratio below y t = 
4 (p t ~ 4 GeV/c) hard-component ratio rAA accurately 
describes evolution of jet structure down to small mo- 
menta (y t ~ 2 or p t R3 0.5 GeV/c) below which system- 
atic uncertainties in soft-component subtraction become 
relatively large. 

The results for low-mass pions (left panel) are consis- 
tent with the hard component in Fig. [28] (right): sup- 
pression at larger y t and compensating enhancement at 
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FIG. 29: Hard component ratios rAA(yt) for five centralities 
of 200 GeV Au-Au collisions (curves) for pions (left) and pro- 
tons (right) [4]. Also shown are p-p data (points, unidentified 
hadrons) compared to hard-component model Hnn- 



smaller y t for more-central Au-Au collisions 4, 45]. The 
solid curve in this panel is derived from the central Au-Au 
points in the previous figure. The points in this panel 
are from the p-p data in Fig. [28^ (left) to provide a refer- 
ence. The result for protons (right panel) is surprising. 
The same factor-five suppression is observed at larger y t 
(near 10 GeV/c), but whereas the energy compensation 
for low-mass pions emerges near 0.5 GeV/c the compen- 
sation for massive protons peaks near y t — 3.5 (jpt ~ 2.5 
GeV/c), substantially above the proton mass (y t ~ 2.7). 

The r AA centrality trend is consistent with that for 
minijet correlations in Ref. I5J: transparency (GLS, 
no "high-p t " suppression) below v — 3 correspond- 
ing to a 50% fractional cross section. The suppres- 
sion/enhancement trend for spectrum hard components 
corresponds quantitatively to the increase of mean frag- 
ment multiplicity inferred from minijet correlations |92j . 

The relation between pion and proton spectrum hard 
components also relates to the "baryon-meson" puz- 
zle [H [55] . The PID hard-component systematics shown 
in Fig. [29] correspond quantitatively to the proton-pion 
puzzle as presented in Ref. [S3] and Fig. [13] [I] . 



Figure 30 (left) shows a comparison between soft com- 
ponent So for 200 GeV p-p and Au-Au spectra (dashed 
curve), a SPS S-S spectrum at 19 GeV (solid points), 
a p-p spectrum at the same energy (open points) and 
a Maxwell-Boltzmann distribution (M-B) with the same 
slope parameter T = 145 MeV [95]. The S-S data have 
been interpreted to indicate radial flow based on devia- 
tion from the MB distribution [9"S] . 

(right) shows radial-flow mean speeds (f3 t ) 
from blast-wave model fits to p-p and Au-Au 



Figure 
inferred 



spectra |97j . Substantial radial flow is inferred for p-p col- 
lisions. It is also notable that radial flow increases most 
rapidly for Au-Au collisions that exhibit transparency to 
low-energy partons (GLS, v < 3) and then less rapidly in 
a more-central interval where jet modification is strong. 

The analysis in Ref. [4] extracted spectrum hard com- 
ponents for identified hadrons (pions and protons) up to 
11 GeV/c from 200 GeV Au-Au collisions. The hard 
components were later described quantitatively with 



FIG. 30: Left: 17 GeV p-p (open points) and 19 GeV fi (solid 
points) m< spectra compared to the soft component Snn of 
the TCM for 200 GeV Au-Au collisions (dashed curve) [95] . 
Right: Radial speed (fit) inferred from blast-wave fits to 62 
and 200 GeV Au-Au collisions [97]. When plotted on par- 
ticipant path length v the relation of the (Jit) trends to the 
minijet sharp transition (hatched band) is interesting. 



a pQCD calculation including modified fragmentation 
functions |45j . The same structure that is interpreted 
as radial flow with a blast-wave spectrum model is de- 
scribed by pQCD as a minimum-bias jet contribution to 
spectra when analyzed with the TCM. The correspon- 
dence is seen in Fig. 30 (right panel). The break in the 
inferred {fit) trend occurs at the sharp transition in mini- 



jet properties (hatched band) reported in Ref. [5]. 



C. Jet-related angular correlations 

The minimum-bias angular correlations from 130 GeV 
Au-Au collisions reported in [6] [33] and shown in Figs. [8] 
and [9] were unexpected and stimulated a followup anal- 
ysis including detailed centrality systematics with high- 
statistics data [5J. 




FIG. 31: 2D angular correlations for peripheral (left, ~ N-N 
collisions) and central (right) 200 GeV Au-Au collisions [Sj- 



Figure 31 shows 2D angular correlations from 85-95% 
and 0-5% central 200 GeV Au-Au collisions [5J. The for- 
mer is approximately equivalent to N-N ss p-p collisions 
and the structure agrees with such measurements [30" | l3"l" ] . 
The data histograms are accurately described by a sim- 
ple model function representing three principal data fea- 
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tures: (a) a same-side (SS) 2D peak, (b) an away-side 
(AS) ID azimuth dipole cos(0a — t) and (c) a nonjet 
azimuth quadrupole cos (2(^»a) [.5- Other model elements 
represent Bose-Einstein correlations, 7-conversion elec- 
tron pairs and projectile-nucleon dissociation. 
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FIG. 33: Left: Nonjet azimuth quadrupole 4q{2D}(6) ob- 
tained with model fits to 2D angular correlations from 62 and 
200 GeV Au-Au collisions 29 . Right: The same data trans- 
formed to conventional measure V2 ■ The curves are generated 
by Eq. (§. 



FIG. 32: 2D model fit parameters for the SS 2D peak ampli- 
tude (left) and AS ID peak amplitude (right) from 200 GeV 
Au-Au collisions [5]. 



Figure 32 shows the amplitudes for same-side 2D (left) 
and away-side ID (right) peak amplitudes [5]. Below 
v = 3 (the sharp transition) the data follow the GLS 
trend expected for transparent Au-Au collisions. The 
scaling with energy in the left panel (factor 1.6) is ex- 
pected for jets based on previous analysis [101] . Note 
that the same jet energy scaling applies to all centrali- 
ties, not just the GLS interval below v = 3. A direct 
mathematical link has been established between these 
jet angular correlations and (pt) fluctuations as shown in 
Fig, p^and Sec.|XEl 



D. Nonjet azimuth quadrupole 



by the relation 
.4 



q = 0.0045R(y?s NN )N bm 



(3) 



where R{^/IWn) = log( A /s^/13.5 GeV)/ log(200/13.5) 
[29] . and e op t is the eccentricity inferred from the optical 
Glauber model. A pQCD prediction based on a color- 
dipole model gives v 2 — 0.02 for p-p collisions, consistent 
with Eq. ^ extrapolated to N-N collisions [55]. 
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Jet-related angular correlations are identified as (a) 
the SS 2D peak and (b) the AS ID dipole. What re- 
mains is (c) the nonjet quadrupole plus some small struc- 
tures (soft component from projectile dissociation, Bose- 
Einstein correlations and conversion-electron pairs) . The 
nonjet (NJ) quadrupole can be identified with "ellip- 
tic flow" measurements [5S]. Jet-related contributions 
to conventional v 2 measurements can be identified with 
"nonflow" 1231. 



Figure 33 (left) shows per-particlc nonjet quadrupole 
amplitude Aq{2D} inferred from 2D model fits to an- 
gular correlations j^Hj- It is notable that the centrality 
trend is approximately a Gaussian on relative impact pa- 
rameter b/bo, where 60 ~ 14-7 fm for Au-Au collisions. 

The right panel shows the same data plotted vs cen- 
trality parameter N part . The Aq data are converted to 
conventional measure v 2 via Aq ~ pov 2 where po is the 
mean single-particle 2D angular density. As noted, half 
of the Au-Au cross section is obscured in the interval be- 
50. The solid and dashed curves are defined 



FIG. 34: Left: SS 2D peak amplitudes from Fig. [32] (left) 
scaled by the number of N-N binary collisions showing strict 
agreement with the expectation for transparent Au-Au colli- 
sions (GLS) over 50% of the t otal cross section [5l 1100] . Right: 
NJ quadrupole data from Fig. [33] (left) replotted on fractional 
cross section. The quadrupole amplitude increases to 60% 
of its maximum over an interval where Au-Au collisions are 
transparent [291H00] . 



low N, 



part 



Figure 34 compares minimum-bias jet data (left panel) 
with the same nonjet quadrupole data plotted on frac- 
tional cross section [5] It is notable that over the 
lower 50% of the fractional cross section where Au-Au 
collisions appear to be transparent to 3 GeV partons 
(minijets follow binary collision scaling or GLS) the non- 
jet quadrupole (right panel) increases to 60% of its max- 
imum value. One can then ask how hydrodynamic flow 
is generated without secondary scattering of partons or 
hadrons [95] 1 100] ? It is also notable that the observed 
energy scaling for jets (1.6, left panel) is similar to that 
for the nonjet quadrupole (1.75, right panel). 
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Figure 35 (left) shows published V2{pt) data for 
minimum-bias 200 GeV Au-Au collisions and identi- 
fied 7r, K and A (points) plotted in a conventional 
format |102H104] . The three curves through data 
are obtained by transformation from a single universal 
quadrupole spectrum model described below and shown 
in Fig. |36| The mass trend for different hadron species 
below 2 GeV/c (so-called mass scaling) is said to confirm 
a hydrodynamic interpretation of v% data as representing 
elliptic flow. 
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FIG. 35: Left: Published PID v 2 (p t ) data for three hadron 
species from Refs. [1U31 1104] plotted in the conventional for- 
mat. The curves are transformed from a common quadrupole 
spectrum described in the text. Right: A data from the left 
panel (open points) are replotted in the form V2(pt)/pt vs 
transverse rapidity y t calculated with the proper hadron mass. 
The solid points are taken from the inset in Fig. [19] [70] . The 
dashed curve is a viscous hydro prediction [105] . 



Figure 35 (right) shows the A data from the left panel 
divided by pt measured in the lab frame and plotted vs 
transverse rapidity y t — ln[(m t +pt)/rrih\ (open points), 
where is the proper hadron mass. The data for all 
three hadrons species in the left panel then have a com- 
mon zero intercept at yt = 0.6 ±0.1. That offset can be 
interpreted as a quadrupole source boost Ay t o (a kind 
of radial flow). The intercept is best defined by more- 
massive hadrons, why the A data are featured. The data 
are consistent with a narrow source-boost distribution 
(an expanding cylindrical shell), whereas (dissipative) 
hydro theory (dashed curve) assumes Hubble expansion 
of a dense bulk medium |105) . The published PID V2(pt) 
data appear to falsify that model. 

The solid points in the right panel are more-recent A v-i 
data for 0-10% central Au-Au collisions shown in Fig. 19 
(inset) [70]. Although the maximum V2 value of the 
more-central data should be smaller, the same zero in- 
tercept near 0.6 is indicated. The correspondence is con- 
sistent with the observation in Ref. |106) that quadrupole 
source boost Ay t o shows no significant centrality depen- 
dence over the most-central 70% of the total cross sec- 
tion for 62 and 200 GeV Au-Au collisions. The more- 
recent negative-going V2 data follow the trend (dash- 
dotted curve) inferred from minimum-bias 2002 and 2004 
RHIC PID data in the quadrupole spectrum analysis of 
Ref. PHI- 



Figure 36 (left) shows the result of a further data 
transformation. «2(Pt) is a ratio - the denomina- 
tor is the single-particle spectrum represented by p{yt) 
that includes a strong jet contribution as demonstrated 
in [3] [55]. The numerator of «2(p*) represents only 
those hadrons that "carry" the nonjet quadrupole. The 
left panel shows the result of multiplying data in the 
V2/pt(lab) format of Fig. 35 (right) by the corresponding 



single-particle hadron spectra in the per-participant form 
(2/N part )p(y t ). The common source boost Ay m = 0.6 is 
again apparent. The solid curves correspond to a com- 
mon source mt spectrum described below. The ratio of 
solid to dashed curves is a simple relativistic kinematic 
factor relating p t (boost) to pt{lab) |102j . 
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FIG. 36: Left: The data from Fig. [35] multiplied by the 
single-particle spectrum on yt for each hadron species reveal- 
ing a common source boost (horizontal offset). Right: The 
data from the left panel transformed to the boost frame (left 
shifted by Ay t o), transformed to m t — and scaled by ex- 
pected hadron abundances (statistical model) to reveal a com- 
mon quadrupole spectrum in the boost frame (solid curve). 



Figure 36 (right) shows the result of four further oper- 
ations: (a) transform from lab to boost frame by shifting 
all spectra in the left panel to the left by Ay t o = 0.6, (b) 
transform spectra in the boost frame from y t to m t , (c) 
multiply spectra by ratio p t (lab) / 'p t (boost) , (d) rescale 
the resulting spectra by the constant factors (1,7,25) in- 
dicated in the panel, consistent with a statistical-model 
prediction for respective hadron species abundances at 
Tchem ~ 150 MeV. Quadrupole data for three hadron 
species are then quantitatively described by a single Levy 
distribution (solid curve) to the uncertainty limits of the 
data. The quadrupole spectrum is cold (T ss 90 MeV) 
and does not correspond to the spectrum representing 
most hadrons (T w 145 MeV). The e+-e" data (dia- 
monds) form a pt spectrum from 91 GeV LEP dijets 
(perpendicular to thrust axis) with a similar low tem- 
perature but smaller exponent n ee « 7.5. 

The data in Fig. [36] represent a minimum-bias av- 
erage over Au-Au centrality. Subsequent analysis has 
determined that quadrupole source boost Ay t o is in- 
dependent of centrality to the uncertainty limits of the 
data [TUB] . We then conclude that from all 200 GeV 
Au-Au PID V2{pt,b) data we obtain two numbers: (a) 
a nonjet quadrupole amplitude Ay t 2 depending on en- 
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ergy and centrality as Eq. ^ and (b) a quadrupole 
source boost A.y t o common to all collision systems. The 
quadrupole spectrum is cold and appears to be univer- 
sal for all Au-Au collision conditions. The solid curve in 
Fig. [36] (right) plus Eq. ^ may then predict all V2{yt, b) 
data for all hadron species. 



E. Fluctuations 

Fluctuations of event-wise mean p t (as a proxy for lo- 
cal temperature), the K/ir ratio (strangeness enhance- 
ment) and other event-wise statistical quantities were 
expected to reveal critical fluctuations associated with a 
QCD phase boundary separating a hadronic phase from 
a partonic (or quark-gluon) phase. A number of fluctu- 
ation measures were proposed |39| . Considerable confu- 
sion emerged as to measure definitions, interpretations 
and the effects of limited detector acceptances. 

a linearized al- 



Figure 
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(left) shows measure Acr Pt: „, 
ternative to r.m.s. fluctuation measure <f> Pt |107) . Further 
study |108j led to the definition of per-particle variance 
difference Aer^ . n |109j directly related to 2D number an- 
gular correlations by an integral equation |110j . 

Figure [37] (left) shows the scale dependence of fluctu- 
ation measure Aap t . n (5r),S4>) with the following proper- 
ties: (a) The bin- wise pt variance at each angular scale 
(St], 5(f>) is compared to a statistical reference, (b) the 
measure is sensitive to p t fluctuations given a bin-wise 
multiplicity n condition (p t : n) 109 . The general trend 
is increase with increasing bin size or scale (5r], Sip). The 
maximum value in this case corresponds to the STAR 
TPC acceptance A77, Acf> — 2, 2ir. Single fluctuation mea- 
surements from detectors having different angular accep- 
tances correspond to different points on that surface. 




FIG. 37: Left: pt fluctuation measure Act? .„ (variance dif- 
ference) plotted vs angle bin sizes (Srj, 5(f)) for 0-5% central 
200 GeV Au-Au collisions |109] . Right: Inversion of the scale 
(bin size) dependence in the left panel to reveal the under- 
lying angular correlations. Fitted AS ID dipole and nonjet 
quadrupole components have been subtracted to emphasize 
the jet-related SS 2D peak with its catenary shape on r/A and 
neighboring negative undershoots. 



Figure 37 



(right) shows the result of inversion of the 
scale dependence of the fluctuation data in the left panel. 
The relation can be described by an integral equation, 



and the integral can be inverted by standard methods to 
reconstruct the underlying angular correlations that gen- 
erated the fluctuations [HOj . We find that (p t ) fluctua- 
tions supposedly related to a QCD phase boundary are 
actually driven by jet correlations. The angular correla- 
tions in the right panel inferred by fluctuation inversion 
were confirmed by direct pair-counting methods. The 
negative undershoots near the SS 2D peak on azimuth 
were also confirmed by direct correlation analysis. The 
result for peripheral Au-Au collisions corresponds to the 
expected jet structure obtained from p-p collisions. Fit- 
ted AS dipole and nonjet quadrupole components have 
been subtracted to emphasize the SS 2D peak as in Fig. [8] 
Reference [101] reports the energy dependence of 
p t fluctuations and angular correlations from SPS to 
RHIC. A simple log(y / s/\T/v/10 GeV) trend is observed 
above 10 GeV, further supporting a jet interpretation. 
Again the similarity with the energy dependence of the 
nonjet quadrupole log(y / ijviv/13.5 GeV) is notable [55] 
(Fig. 34 1 . The two intercept points on energy are con- 
sistent within systematic uncertainties and are kinemat- 
ically consistent with the first significant penetration of 
projectile nucleons significantly below valence quarks on 
momentum fraction x. 



XI. DISCUSSION 

Two distinct approaches have been applied to RHIC 
heavy ion data, (a) Evidence for a conjectured quark 
gluon plasma has been sought in the form of "signals" re- 
covered by model-dependent analysis methods designed 
to detect them. Priority has been given to results that 
seem to confirm QGP formation, (b) Significant struc- 
tures in yields, spectra and correlations are revealed by 
differential methods not depending on a priori physical 
models. Significant structures are represented by simple 
functional forms, and the systematic variation of model- 
function parameters with collision conditions is obtained. 
Model-parameter systematics are then compared with 
physical models, giving first priority to QCD descrip- 
tions of elementary collisions. The two approaches have 
led to quite different interpretations of RHIC data. In 
this section we consider the differences and their possible 
implications. 



A. Initial-state energy density 

Estimating the initial conditions (mid-rapidity gluon 
and energy densities) for A-A collisions at RHIC is said 
to be essential to support claims of densities and tem- 
peratures sufficiently high to produce deconfined quarks 
and gluons according to LQCD and to support a hydro- 
dynamic description in terms of large initial pressure gra- 
dients 121] . The initial-state energy density is convention- 
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ally estimated by the Bjorken formula 



(4) 



1 dE t 
tA drj 

1 , x dn ch 

W T(m t ) — ; 

tA x ' drj 
ex ^[1 + 0.1(1/ -1)], 

where the last line assumes overlap area A oc Np^ t 
and we invoke the two-component expression for dn c h/dr] 
from Ref. [8]. Reference [43] estimates (e) sa 15 GeV/fm 3 
for central 200 GeV Au-Au collisions based on (mt) ~ 
0.85 GeV/c 2 , seemingly well above the threshold for 
producing a QGP. But what is the energy density at 
the sharp transition reported in Ref. j5]? In that case 
Npart ~ 50 rather than 380 and v 3 rather than 6. 
Based on the centrality scaling above the energy density 
is reduced by only a factor 0.4, to 6 GeV/fm 3 , nominally 
well above the critical density required for QGP forma- 
tion. But spectrum and angular correlation data indicate 
that Au-Au collisions are still transparent at that cen- 
trality, with no evidence for parton or hadron secondary 
scattering [5]. Applying the Bjorken formula assuming 
that all produced E t is subject to thermalization may 
greatly overestimate the relevant energy density. Much 
of the final-state E t may be carried by hadrons (and par- 
ent partons) that experienced no secondary interactions. 



B. Initial conditions, the TCM and minijets 

The HEP two-component model based on projectile 
dissociation and minijet (minimum-bias jet) production 
has been discounted as a description of initial conditions 
at RHIC based on comparisons between HIJING simu- 
lations and measured hadron production systematics at 
mid-rapidity [3 (2TJ [47] . Reference [11] states that "...the 
available RHIC experimental data... can already rule out 
the simple two-component model without nuclear modifi- 
cation of the parton distributions in nuclei...." The basis 
for that conclusion is a TCM for p-p collisions given by 



dn r 



dr) 



n h - 



a jet 
0~inel 



(5) 



where n s = 1.6 and — 2.2 (dijet fragment multiplic- 
ity) are assumed [11]. The measured charge density for 
200 GeV NSD p-p collisions is about 2.5, implying that 
the cross-section ratio must be about 0.4. Given that 
o~inei ~ 60 mb or more in that reference the calculated 
a j e t must be at least 25 mb, more than half the actual 
p-p inelastic cross section 42 mb. It is not clear how a 
pQCD calculation can determine such a nonperturbative 
cross section. 

The generalization of Eq. ^ to A-A collisions is given 
by ® 



2 dn r 



where n pp ~ 2.5 represents the 200 GeV p-p NSD r\ den- 
sity and x = nhO-jet/n-ppO~inei is the same parameter ap- 
pearing in the TCM from Ref. 0. The value x = 0.35 
inferred from Ref. |11| can be compared with x i=a 0.1 
inferred empirically from more-central 200 GeV Au-Au 
data. Reference [llj concludes that minijet production 
must be suppressed by "shadowing" in Au-Au collisions. 

But the HIJING (PYTHIA) Monte Carlo characteri- 
zation of p-p collisions is quantitatively inconsistent with 
p-p data. Reference [IB] reported a TCM for p-p spec- 
tra from which a dijet production probability could be 
inferred corresponding to a 2.5 mb dijet total cross sec- 
tion [IS] (Fig. 27 of this paper). That cross section is 
consistent with a parton spectrum lower bound near 3 
GeV (dijet energy Q — 6 GeV) and also consistent with 
the original minijet observation by UA1 [llll and theo- 
retical descriptions thereof [1121 1113| . In Eq. Q the cor- 
responding value is x ~ 0.03 (GLS extrapolation), more 
than ten times smaller than assumed in Ref. [11) . The jet 
cross section is overestimated in part because the steep 
power-law parton spectrum is assumed to extend down 
to po — 2 GeV, well below 3 GeV inferred from spectra 
and correlations [!45, 46| [92] . The fragment production in 
more-central Au-Au collisions is actually underestimated 
by extrapolation from measured p-p data (GLS). Modi- 
fied fragmentation leads to the larger x 0.1 in more- 
central collisions (Figs. 26 and 28 of this paper). The 



parton spectrum cutoff, and therefore jet cross section, 
is nonperturbative and must be inferred from data. 



N, 



part 



drj 



,[l + x(v-l)], 



(6) 



Figure 11 (left) compares HIJING predictions (dotted 
lines) with data (points) at 20 and 200 GeV. The HI- 
JING slope at 200 GeV (hard component) is inconsis- 
tent with the data, one basis for discounting the mini- 
jet TCM. But setting aside the Monte Carlo simulations 
the hard component inferred from p-p data |46j is ten 
times smaller and underpredicts the Au-Au data. Within 
the lower 50% of the cross section (below N part ~ 50, 
not shown in those plots) the p-p extrapolation accu- 
rately describes more-peripheral Au-Au data [52] • Above 
that point (sharp transition [5 ) fragment production in- 
creases by a factor 3-4 as jet production is modified by 
the medium |3||35]. The two-component minijet model 
accurately describes hadron production (yields, spectra 
and correlations) in terms of pQCD parton scattering and 
nonperturbative (measured) fragmentation to jets. The 
relation between HIJING simulated yields and correla- 
tions and Au-Au measured correlation data is discussed 
in detail in Sec. VIII-I of Ref. ED . 



C. The CGC vs HIJING and data 

Reference [5T] concludes "If the initial conditions were 
not well constrained at RHIC, then the conclusion that 
a QGP was formed could not be sustained. [...] ...the 
surprising very weak centrality and beam energy depen- 
dence... is most satisfactorily explained and predicted by 
the CGC... This is one of the strongest lines of empirical 
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evidence... that the CGC initial state... is formed...." The 
reference is to HIJING simulations compared to Au-Au 
data as in Fig. [IT] of this paper. As noted, HIJING with 
parton spectrum lower bound po — 2 GeV greatly over- 
predicts the gluon density or minijet production com- 
pared to what is inferred directly from p-p collisions. A 
CGC estimate based on a saturation scale Q s « 1-1.5 
GeV is even larger, leading to inference of an opaque (glu- 
onic) medium. And the centrality variation of the CGC 
estimate is seen to disagree strongly with the measured 
centrality trend when observed over the entire centrality 
interval as in Fig. 26 not just more-central collisions [55] , 



D. Jet contributions to spectra and correlations 

Blast-wave (BW) fits to pt spectra are employed to 
infer bulk-matter properties such as temperatures and 
radial flow speeds (Figs. [2] and 12). The assumption 
that spectrum structure actually represents such quanti- 
ties can be questioned [SS]. In Fig. 28 (right) we observe 
that a pQCD description of the pion spectrum hard com- 
ponent provides an accurate and quantitative description 
for all Au-Au centralities. In Fig. |30] (right panel) we find 
that the radial speed inferred from BW fits is rising most 
rapidly over a centrality interval (GLS) where Au-Au col- 
lisions appear to be transparent according to jet angular 
correlations [5]. The nonzero (f3t) inferred for 200 GeV 
p-p collisions characterizes a pt spectrum shape similar to 
that for 19 GeV fi collisions. Does the departure from a 
M-B distribution signal radial flow or a universal aspect 
of longitudinal nucleon or parton fragmentation? 

Baryon-meson spectrum ratios (Fig. 13 ) are inter- 



preted to favor hadronization from a thermalized par- 
tonic medium by recombination of constituent quarks. 
However, the spectrum hard components in Fig. [29] pro- 
duce the same B /M anomaly but can be associated with 
parton fragmentation to jets [4]. The most-central pion 
data from Fig. [29] are the same data (points) described 
by pQCD (solid curve) in Fig. 28 (right). 



Jet quenching is inferred from spectrum ratios (Fig. 14 ) 
that reveal jet structure only above 4 GeV/c. In con- 



trast, spectrum hard-component ratios (Fig. 29 ) reveal 
jet structure down to 0.5 GeV/c {y t ~ 2). Hadron sup- 
pression at 10 GeV/c is accompanied by much greater 
enhancement at 0.5 GeV/c (for pions) consistent with a 
simple mechanism for modified fragmentation 45 . The 
low-p t enhancement remains within the resolved jet as in 



Fig. 31 5, 92, 95J, with jet parameters as in Fig. 32 



Jet correlation structure on ID azimuth is inferred 
by ZYAM subtraction of a combinatorial background 
with certain trigger-associated pt cuts imposed (Figs. 22 
and 23). The ZYAM subtraction method assumes that 



SS and AS jet peaks do not overlap and that published v 2 
data are unbiased. However, we observe that minimum- 
bias SS and AS jet peaks overlap strongly on azimuth 
and that much of the published i>2 data include strong 
biases from jet contributions (nonflow) [2H1 1106] leading 



to (a) oversubtraction of the constant offset and sinusoid 
and (b) artificial suppression and distortion of ZYAM- 
inferred jet structure, including an AS double peak and 
inference of "Mach cones" |76j . 

Jet correlation structure on 2D (rj, <f>) with trigger- 
associated p t cuts imposed reveals substantial changes 
in the SS 2D jet peak shape with increasing Au-Au cen- 
trality (Figs. 24 and 25). The ^-elongated SS peak in 



more-central collisions is separated into a "jet-like" part 
narrow on An and a "ridge-like" part broad on An. It 
is notable that both parts share the same azimuth width 
and the same energy systematics. But the "ridge" is 
attributed to nonjet mechanisms. Minimum-bias 2D an- 
gular correlations (no trigger-associated p t cuts) show a 
single monolithic SS peak well-described by a 2D Gaus- 
sian for all centralities [3] . The peak volume corresponds 
to pQCD predictions for jet systematics [92!. p t angu- 



lar correlations from fluctuations (Fig. 37 1 indicate jet 



like elevated pt values in the elongated SS peak region 
on r]A smoothly connected to the central region with no 
discontinuity. The SS 2D peak structure resulting from 
trigger-associated p t cuts represents a small fraction of 
the minimum-bias jet peak. 



E. The nonjet quadrupole and elliptic flow 

"Elliptic flow" V2 is conventionally measured by several 
nongraphical numerical methods (NGNM) as in Fig. 15 
The range of results for different analysis methods is the 
basis for estimating systematic uncertainties. The source 
of such variation may be dominated by jets [106) . An 
alternative method based on model fits to 2D angular 



correlations as in Fig. 33 returns reproducible data that 
distinguish accurately between jet structure and a nonjet 
azimuth quadrupole |29j . 

So-called mass scaling of v 2 {pt) below 2 GeV/c is said 
to confirm a hydrodynamic phenomenon. However, in 



Fig. 34 the nonjet quadrupole amplitude is observed to 
rise to 60% of its maximum value over the GLS interval 
(lower 50% of total cross section) where Au-Au collisions 
are transparent to 3 GeV partons: the SS peak amplitude 
scales with the number of N-N binary collisions while the 
peak shape is unchanged from p-p collisions |100) . 

In Figs. [17] through [19] attempts are made to demon- 
strate so-called constituent-quark scaling above 2 GeV/c 
interpreted to demonstrate flow of a quark-gluon medium 
in which constituent quarks play a prominent role. But 
Figs. 35 and 36 demonstrate that V^ipt) f° r au hadron 



species can be predicted from a single universal spec- 
trum shape with these characteristics: (a) fixed source 
boost Ay m w 0.6, (b) temperature T 2 « 90 MeV, (c) 
Levy shape with exponent n s» 14.5. Similar analy- 
sis of centrality dependence reveals that Ay t o does not 
vary significantly with Au-Au centrality over the interval 
70-0% [106]. In Fig. [35] (right) we find that the boost 
distribution inferred from VzijPt) data is very different 
from that assumed in hydro predictions for a Hubble- 
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expanding medium. 

So-called V2 "scaling" methods tend to obscure the 
sought-after phenomenon. Hydro theory can only pre- 
dict the velocity (rapidity boost) profile of a flowing 
medium. Hadron details must be supplied by a sep- 
arate model. If i>2 data can be manipulated to reveal 
the hadron source boost then direct comparisons can be 



made between data and hydro. In Figs. 35 (right) and 



[36] the common source boost for all hadrons can be read 
directly as intercept Ay t o with good sensitivity. On p t 
the intercept points for different hadron masses become 
Pto = rrih sinh(Ayto) ~ m^Ay t o ps 0.6m/, (so-called mass 
scaling). If the data are transformed to "kinetic en- 
ergy" KE t (or equivalently transverse mass) the inter- 
cept points become m t0 — = m/ l [cosh(Ay t0 ) — 1] ~ 
0.5m/j(Ay t o) 2 ~ 0.18m^. The sensitivity to boost is fur- 
ther reduced by a factor 3. Plotting on y t maximizes 
sensitivity to hydro boost predictions while plotting on 
ffh ~ m h minimizes any sensitivity to the nominal theory 
Certain differences between V2 methods interpreted to 
indicate nonflow and/or flow fluctuations as in Fig. [l6] 
(right) are quantitatively predicted by the properties of 
the SS 2D peak from jet correlations The i>3 (tri- 

angular flow) data in Fig. 21 correspond exactly to the 



SS 2D jet peak in 2D angular correlations as in Fig. 31 
A similar correspondence is shown for other "higher har- 
monics" |551ITT3]. 

V2{pt) is a ratio — the denominator is the single-particle 
spectrum p(pt) including a strong jet contribution over 
a broad p t interval; the numerator includes as a factor 
the quadrupole spectrum P2(pt) with no significant jet 
contribution |102j . From Fig. 36 we learn that pzipt) 



is a cold spectrum very different from the single-particle 
spectrum describing most hadrons. It is possible that 
the nonjet quadrupole is actually "carried" by a small 
fraction of all hadrons. Ideal hydro calculations implicitly 
assume a common spectrum for almost all hadrons which 
should cancel in the ratio. The falloff from ideal hydro 
in Fig. [18] and similar results reflects the differences in 
the two spectra combined in ratio. Only the numerator 
of v-zipt) has significance for "elliptic flow." Hydro theory 
cannot in principle describe the V2 ratio over significant 
Pt intervals. 



F. Gluons in dijet and quadrupole production 

The possible relations among soft gluons released dur- 
ing projectile nucleon dissociation, midrapidity dijets and 
the nonjet quadrupole are intriguing. If one sets aside 
a hydrodynamic mechanism the sharply-peaked rj struc- 
ture in Fig. 4] (right) suggests a low-x glue origin for the 
nonjet quadrupole [35]. In Fig. 20 (lower right) the r] 
dependence is shown to scale closely with beam rapidity. 

The energy dependence of dijet production is deter- 
mined largely by the depth on x of the nucleon PDF 
probed by a given collision energy. Once the longitudi- 
nal gluon flux is so determined dijet production by scat- 



tering is determined by the dijet energy scale Q with 
fixed lower bound near 6 GeV. Soft and hard hadron 
production should then vary similarly with energy. We 
observe that the nonjet quadrupole scales with energy in 
essentially the same way, suggesting low- a; glue as a com- 
mon origin for soft, hard and quadrupole. pQCD calcu- 
lations based on a color-dipole model predict p t -integral 
V2 = 0.02 for p-p |99j consistent with peripheral Au-Au 
collisions [33]. Thus, dijets and nonjet quadrupole may 
have a common gluonic origin but represent two limit- 
ing cases of QCD: short- wavelength and long- wavelength 
gluonic radiation [2"2"] . 



XII. SUMMARY 

The RHIC accelerator complex has been a major tech- 
nical success. The development of computing, software 
and analysis methods to deal with petabytes of particle 
data have required an enormous effort by well over 1000 
persons. The product has been a massive inventory of 
analysis results reported in hundreds of published papers 
relating to a very complex question: what physical mech- 
anisms dominate ultrarelativistic heavy ion collisions and 
what state of matter (if any) results? 

Strong theoretical arguments have been advanced to 
support the claim that a strongly- coupled quark gluon 
plasma (sQGP) is formed in more-central Au-Au colli- 
sions with a remarkably small shear viscosity, warranting 
the designation "perfect liquid." The principal experi- 
mental observations invoked are large-amplitude elliptic 
flow v 2 data described to good approximation by ideal hy- 
drodynamics (indicating small dissipation) and strong jet 
quenching, apparent five-fold reduction of jet fragment 
yields by a dense QCD medium nearly opaque to colored 
quarks and gluons. Other evidence includes strong radial 
flow inferred from p t spectra and mid-rapidity hadron 
yields indicating large initial gluon (color glass conden- 
sate) and transverse-energy E t densities that may drive 
the hydrodynamic flows. 

However, a two-component model (TCM) of A-A col- 
lisions based on measured properties of p-p collisions is 
found to describe accurately the centrality trends from 
Au-Au collisions across the more-peripheral half of the 
total cross section. And the systematics of jet-related 
structure in the more-central half is still described in part 
by perturbative QCD (pQCD). Within the same more- 
peripheral interval a "quadrupole" measure, alternative 
to V2 but describing the same correlation structure iden- 
tified as elliptic flow, increases to 60% of its maximum 
value within peripheral Au-Au collisions where no sec- 
ondary scattering of partons or hadrons is observed. And 
in more-central collisions the same simple quadrupole 
trend remains unchanged, although jet structure (parton 
fragmentation) is strongly modified there. 

Differential analysis of ^(pt) data reveals a quadrupole 
source boost distribution and unique quadrupole m t spec- 
trum that may be compared directly with hydro predic- 
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tions. The comparison contradicts hydro models based 
on Hubble expansion of a bulk medium as assumed for 
sQGP. The quadrupole spectrum is quite different from 
the spectrum describing most final-state hadrons, sug- 
gesting that the quadrupole identified with elliptic flow 
is actually carried by a small minority of all final-state 
hadrons. Other results for hadron production, p t spec- 
trum structure and jet correlations also conflict with the 
sQGP paradigm. 

The correct model for high-energy heavy ion collisions 
may lie on a continuum between production of projec- 
tile and scattered-parton fragments (hadrons) described 
by the two-component model and late hadronization of a 
thermalized weakly-coupled QGP as limiting cases. The 
sQGP conjecture already retreats along that continuum 
from the asymptotic weakly-coupled QGP. However, fur- 
ther paradigm change may be necessary. Three major 



questions emerge from a review of RHIC experimental re- 
sults: (a) What is the mechanism for the nonjet azimuth 
quadrupole? (b) What mechanism modifies parton frag- 
mentation to jets in more-central Au-Au collisions? (c) 
Is there any necessary role for hydrodynamics in RHIC 
collisions? Some of the exceptions to the standard heavy- 
ion paradigm emerging in recent years suggest that QCD 
is a richer field theory than was previously imagined. 
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